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DECLARATION OF SHIGEO T. OYAMA ("TED OYAMA") 
UNDER 37 C.F.R. §1.132 

Shigeo T. Oyama declares as follows: 

1 . I am the inventor of the above-identified application. I am also an expert in the 
field of catalysts qualified to provide experimental evidence and opinion evidence 
directed to such matters at the level of skill of one of ordinary skill in the art and 
obviousness to one of ordinary skill in the art. I have attached hereto a one page, 
short form resume, to establish my credentials as an expert. As demonstrated by 
the resume, I hold the degree of Ph.D. in Chemical Engineering; I have been 
involved in research in Chemical Engineering, Chemistry, an Catalysts for over 
two decades; I am a titled professor at Virginia Polytechnic Institute & State 
University; I am a noted author on over 140 refereed publications; I have edited 
several books; and I am on the Editorial Board of the Journal of Natural Gas 
Chemistry. 



2. I have read and understand the application as originally filed, the amendment to 
be concurrently submitted with this declaration, the office action mailed July 30, 
2004, and each of the references cited by the Examiner. It is my opinion that the 
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claimed invention, particularly as amended in claims 1 and 7, would not be 
anticipated by or be obvious to one of ordinary skill in the art over U.S. Patents 
4,454,246 and 4,359,406 to Fung, the Journal of Catalysis article (XP009009979), 
or U.S. Patent 4,367,137 to Antos. 

A. Claimed invention 

The claimed invention can best be understood with reference to 
independent claims 1 and 7, as amended, which recite as follows: 

1. (Currently Amended) A catalyst comprising: 

a metal phosphide complex having the formula MP,, wherein M is selected 
from the group consisting of V, Cr, Mn, Fe, Co, Ni, Nb, Mo, Ta, and W, and 
wherein x ranges from about 0.1 to about 10; and 

a high surface area support of at least 50m 2 /g . wherein the metal phosphide 
complex is dispersed on the high surface area suppo rt, wherein said hi> h surface 
area support is selected from the g roup consisting of carbon, silica, titania. thnria , 
magnesia, zirconia. kaolin, benton ite. kieselg nhr . zeolites, and combinations 
thereof 

7. (Currently Amended) A metal phosphide catalyst comprising: 

a metal phosphide complex having the formula A a B b P y , wherein A and B 
are each selected from the group consisting of V, Cr, Mn, Fe, Co, Ni, Nb, Mo, Ta, 
and W, wherein the sum of a and b is 1, the ratio of a and b ranges from about 0.01 
to about 100, and y ranges from about 0.1 to about 10; and 

a high surface area support of at least 50m 2 /g wherein the metal phosphide 
complex is dispersed on the high surface area suppo rt, wherein said hi P h surface 
area support is selected from the gr oup consisting of carbon, silica, titania. thoria , 
magnesia, zirconia, kaolin, bentonite. kieselguhr. zeolites, and combinations 
thereof. 



is noted that the claims specify particular metals ( V. Cr. Mn. Fe. Co. Ni. 
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Nb, Mo, Ta, and W) which are "base" metals. In contrast the Fung and Antos 
patents are directed to "noble" metal phosphides. As will be discussed below in 
more detail, neither Fung nor Antos describe procedures which would yield metal 
phosphides from the specified metals (V, Cr, Mn, Fe, Co, Ni, Nb, Mo, Ta, and W). 
In view of this, these references do not anticipate and would not make obvious to 
one of ordinary skill in the art that metal phosphide catalysts of the particular 
metals (V, Cr, Mn. Fe. Co. Ni. Nb. Mo. Ta. and WV A person of ordinary skill in 
the art to which the present invention pertains would typically hold a Ph.D. and 
typically with 5-15 years research experience after graduation. He or she would 
have access to a variety of trade journals and reference books in the fields of 
chemistry, chemical engineering and catalysis, and would be skilled at reviewing 
the data presented in those references. 

Furthermore, the claims specify particular high surface area supports which 
include carbon, silica, titania. t horia. magnesia, zirconia. kaolin, hentonite , 
kieselguhr. zeolites, and combinations thereof . Unlike claims 2 and 8 of the 
application, which have now been canceled, claims 1 and 7, as amended, do not 
specify alumina as the support. Alumina is the support used in the Nozaki 
(Journal of Catalysis 79, 207 (1983)) reference identified as (XP009009979). As 
will be discussed in more detail below, alumina is a poor support for phosphides 
because of its strong interaction with phosphorous to form aluminum phosphate 
(A1P0 4 ) and necessitates the use of higher temperatures for reduction. 

B. U.S. Patents 4,454,246 and 4,359,406 to Fung 

These patents describe the preparation of metal-phosphorous compounds 
deposited on a high surface area support. While the patents discuss Group VIII 
metals generally, and Group VIII noble metals specifically, it is noted that the 
synthesis method described is inappropriate for the non-noble metals. 
Specifically, the reduction temperature employed in Fung is low (300-500°C), 
which is appropriate for noble metals, but not for base metals (e.g., V, Cr, Mn, Fe, 
Co, Ni, Nb, Mo, Ta, and W). 
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In Fung, there is described a procedure for impregnation of a support with 
soluble metal compounds, a phosphorous source, and reduction of the compounds 
directly in hydrogen, without high temperature calcination. Note particularly 
Examples 1-2 of the Fung references where there is described low temperature 
drying at 120°C before reduction. Although, this method works for the 
preparation of PtP 2 because Pt can be reduced easily. It does not work for non- 
noble metal phosphides. 

In sharp contrast, in the present invention, the metal precursors are first 
calcined at high temperature to ensure that the components are combined in the 
form of a phosphate, and then the reducing is performed. If this is not done, the 
components are not combined , and metal phosphide catalysts will not be 
produced. To demonstrate this, the following experiments was performed: 

Method of Reduction of Samples 

Phosphide samples supported on a silica support were prepared by 
impregnating the support with aqueous metal nitrate and ammonium phosphate 
solutions by the method of incipient wetness impregnation. Prior to use, the silica 
was dried at 393 K for 3 h and calcined at 773 K for 6 h. The incipient wetness 
point of the silica was found to be 2.2 cm 3 g'. After impregnation, the powders 
were dried at 393 K for 3 h and either used as such or calcined at 773 K for 6 h. 
The reduction of the samples was carried out using temperature-programmed 
reduction (TPR). The reduction was carried out in a U-shaped quartz reactor 
placed in a furnace controlled by a temperature programmer (Omega Model CN 
2000). A portion of the exit gas flow was sampled through a leak valve into a 
mass spectrometer (Ametek/Dycor Model MA 100) for following the progress of 
the reaction. The temperature was monitored by a local chromel-alumel 
thermocouple placed in a thermowell near the center of the reactor bed. The H 2 
flow rate was set at 1000 mol s° (1500 cm 3 min 1 ) per gram sample. The heating 
rate was 1 Kmin" 1 . 
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Results 

Reduction of Components 



Compound 


Reduction 
temperature 
separate 
components 


Result 


Reduction 
temperature 
combined 
components 


Result 


FeP/Si0 2 


570-700 K 
(300-430°C) 


Metal 


930 K 
(660°C) 


phosphide 


CoP/Si0 2 


540-590 K 
(270-320°C) 


Metal 


830 K 
(560°C) 


phosphide 


Ni 2 P/Si0 2 


530-650 K 
(250-280°C) 


Metal 


820 K 
(550°C) 


phosphide 



As can be seen from the above table, the precursors employed are metal salts like 
nitrates, and a phosphorous source like ammonium phosphate. The metal has a 
positive charge which is balanced by the negative nitrate ion and the phosphate 
has a negative charge which is balanced by the positive ammonium ion. Because 
the charges are balanced, there is no need for the metal to combined with the 
phosphate (so the result is a "metal"). However, in the calcination step, the nitrate 
ion and the ammonium ion are destroyed, to form gaseous nitrogen oxides and 
ammonia which are released from the solid. This leaves behind the positive metal 
ion and the negative phosphate species which combine in order to balance their 
charges. 

The table above demonstrates that direct reduction of components, as 
contemplated by the Fung references, does not work for the base metals recited in 
the claims. Because the metal component is not combined with phosphorous, the 
metal is reduced to a metallic state at a low temperature, and then sinters (grows 
particle size) as the temperature is raised, precluding its combination with 
phosphorous. 

One of ordinary skill in the art would recognize that the results in the Fung 
patents are limited to Pt phosphides (this being the only metal where test results 
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are presented). Further, one of ordinary skill in the art would recognize that Pt is 
a noble metal, and is chemically different from the base metals recited in the 
claims. Thus, one of ordinary skill in the art would recognize that Fung's data 
could not be extended to all group VIII metals. This is verified by the results in 
the Table above. 

C. U.S. Patent 4,367,137 to Antos et al. 

i 

This patent describes an acidic catalytic composite composed of a carrier 
metal supporting a platinum group metal, cobalt, tin, phosphorous, and halogen. It 
is noted that nowhere in Antos is it demonstrated that a phosphide is formed. 
Rather, only a brief statement in the text alludes to this, but acknowledges 
uncertainty. See particularly column 14 where it is stated: 

"Despite the fact the precise form of the chemistry of association of the 
phosphorous component with the catalytic composite of the present invention is 
not known , it is believed to be most intimately related with the cobalt component 
of the catalytic composite, although it may be associated physically and/or 
chemically with the carrier marerial and/or the platinum group and tin 
components"... "it is believed that best results are obtained when the phosphorous 
component exists in the catalytic composite substantially in the form of a 
phosphide with the cobalt component." 

In fact, the procedures used by Antos would not yield a metal phosphide as 
claimed by the present application. This is because the temperature used in 
Example II of Antos (reduction of the Co and P containing composite at 975 °F 
(524°C), is too low for phosphide formation. As evidence of this, attached hereto 
are the following: 

1) Goodenough, J. Appl. Phys. Vol. 40, p. 1250 (1969). This document reports 
that phosphide formation takes place been 850-900°C between metal and 
phosphorous. See also, Schonbert, Acta Chemica Scandinavica 8:226-239 (1954) 

2) Wang et al. /. Catalysis 208-321-331 (2002). This document, on which I am 
noted as joint-author, shows that with Co/P, major reduction took place at 820 K 
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(550°C), but completion required 690°C. (See Figure 4 on page 325). This 
corresponds to Example 5 of the present application where 690°C was used. 

Based on these reported results, as well as the lack of data and analysis in 
Antos, and Antos' own statements of uncertainty, it should be concluded that 
Antos does not suggest a methodology which would yield metal phosphide 
catalysts which fall within the ambit of the claims of the present application. 
Thus, one of ordinary skill in the art would not find Antos as suggesting the 
claimed invention. 



D. Nozaki et al, Journal of Catalysis 79, 207 (1983) 

In this paper, the authors report catalysts composed of Ni 2 P, Co 2 P, and FeP 
deposited on an alumina support. The catalysts were obtained from phosphate 
precursors by reduction in hydrogen. 

One of ordinary skill in the art will recognize that the use of other supports, 
as set forth in independent claims 1 and 7, provide certain new and unobvious 
advantages over alumina supports. In fact, alumina is a poor support for 
phosphides, as it interacts very strongly with the phosphorus forming an aluminum 
phosphate (A1P0 4 ) and necessitates that use of high temperatures nf reduction 
The attached figures, which are taken from referenced publications on which I am 
noted as a joint author, show reduction provides for the different catalysts obtained 
by following the production of water as a function of temperature. In all cases, the 
precursor is a phosphate which undergoes reduction to a phosphide. 

Attached Figure 1 shows the reduction profile of pure MoP and MoP/Si0 2 
(see, Clark et al., J. Catal. 207,256, (2002) attached). Figure 2 shows the 
reduction provide for MoP/Al 2 0 3 (see, Clark et al., J. Catal. 218, 78 (2003) 
attached). As can be seen, the reduction for MoP/Si0 2 requires only 820 K 
(547°C), while that for MoP/Al 2 0 3 necessitates 1 170 k (900°C). Published work 
also shows that this is because of the formation of the A1P0 4 phase as a result of 
the strong interaction between Al and P. 

Similar results are found for Ni 2 P on other supports claimed in the present 
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application. Attached Figure 3 compares the synthesis ofNijP supported on two 
zeolites (KUSY. potassium promoted ultra-stable Y) and MCM-41 , to Ni 2 P 
supported on silica and Ni 2 P supported on Al 2 0 3 . The reduction teimperaiures for 
the KUSY, MCM-41, and SiOj supports are between 860*880 K (590-610° C), 
-whereas mat for Al 2 0 9 support is at 1270 K (1000 9 C). The presence of Ni^P has 
been confirmed on all these supports by x-ray diffraction and x-ray absorption 
measurements, this work shows that supports other than alumina can produce the 
phosphide at much lower temperatures. This is a non-obvious, and highly 
important finding. 

Attached Figure 4 shows die results of preparing MoP, WP, and Ni a P using 
a carbon support. Again, the reduction temperatures are low, ranging from 810- 
860K(540-590°C)- 

3. 1 hereby declare that all statements made herein of ray own knowledge are true 
and that all statements made on information and belief are believed to be true; and 
further that these statements were made with the knowledge mat willful false 
statements and the like so made are punishable by fine or imprisonment, or both, 
under 18 U.S.C. 1001 and that such willful false statements may jeopardize the 
validity of the application or any patent issued thereon. 
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JOURNAL OF APPLIED PHYSICS 

Transition Metal Compounds 

J. B. Goodenough, Chairman 

CrystaUographic and Magnetic Properties of Solid Solutions of the 
Lty Phosphides MlP, M=Cr, Mn, Fe, Co, and Ni 

R. Fruchart, A. Roger, and J. P. Senateur 
C.E.C.M., 15 Rue G. Urbain, 94 Vitry sur Seine, France 

Several phosphides with the chemical formula ™<g$^£3££ 
and Ni,P. The structure of Co,P is orthorhombic. These two struct ure i are reUt edm an 
£y subce.l consisting of a tetrahedral-site and gg^^Jg^ Kh M"'^ "hibit 
that solid solutions between two he«gonal end n ^ b ^" C ^ w n„ m p, ete solid solubUity with NW>. 
orthorhombic structures at intermediate f ™ ^ number of 3d electrons, 
Curiously, h«agonal and orthorhoinbic sy^ ^ 
Ni ! P(he X )-Co,P(ortho)-Fe J P(h«)-FeMnI ortho) Mn^lhex^ P Cr substituting 

composition and ^Messbauer studies preferentially for tetrahedral- 
for pyramidal-site iron m FeMnP and tecrr wm " 1> * _ f errom aenetic, intermediate phases 
site konTn FeCoP and FeNiP. While neither Co,P ^J^^^SsSrC and 3.03 „ B /mole- 
are, the Curie temperature and magneUzat.on ^reachmg ^the ^ ( m„ 0 . 7 Co,.,),P. The 

cule for MnCoP. Metamagnetism k'*^* SuriTSSL* near FeCoP. This system 

Co.P-Fe.P system is also ferromagnetic wuh » ^"^^SL^ in Curie temperature, across 

atom ratio near that for the appearance of metamagnetism. 
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I. INTRODUCTION 

Several phosphides of the first-row transition metals 
are ferromagnetic: for example, Fe,P >■* Fe 2 P,i ' < and 
MnP s Near 50°K, MnP exhibits an antiferromagnetic 
to ferromagnetic transition,^" the antiferromagnetic 
phase consisting of two parallel helices." - Other phos- 
phides, such as Mn,P » and Mn 2 P,« are antiferromag- 

^Studies of MnP-FeP solid solutions show that FeP 
is not ferromagnetic even at very low (1.5 K) temper- 
atures' 41 ^ 18 MnP-MnAs soUd solutions are mag- 
netically complex, exhibiting several magnetic tran- 

. H . Le Chatelier and S. Wologdine, Compt. Rend. 149, 709 

Haughton, J. Iron Steel Inst 115 417 (1927) . 
tihn Chiba T Phys. Soc. Japan 15, 581 (I960). 
4 IT J P Meyer and M. C. 6 adeville, J. Phys. Soc. Japan 17, 

^iffiweeny and A. B. Scoot J. Chem. Phys », 917 (1954) 
e B F. Stein and R. H. Walmsley, Phys Rev^. 148, 933 (1966) . 
t G P. Felchert, J. AppL Phys 37, 1056 (lW. 
• J Komatsubara, K. Kinoshita, and E. Hirahara, J. rnys. 

^A^^Sart, Compt. Rend. 264, 508 (1967, 

'Hubert Jr. and D.H.Ridgley, Phys. Rev. 135, 1099 . 2S1 fl 967) • , ... . ,y 



sitions Substitutions of boron for phosphorous in J 
F e p24.» have shown that the metaUoid electrons are . 
donors to the 3d bands of iron, whereas substitutions of. ' 
chromium, 2 ' nickel or cobalt" " give changes of mag^ 
netization with composition similar to those found * 
for the borides Af 2 B and MB. j& 
This work is concerned with the magnetic propertied 
of solid solutions of phosphides of the type JfiP| 
Nowotny and Henglein" have presented a brief struc. 
tural survey of a few solid solutions among MiP pnosf 
phides. Not all the phosphides 3f 2 P have &e saine, 
crystal structure*^; some, like Mn 2 P and Ni 2 P, are 

»J. B. Goodenough, D. H. Ridgley, and W. A. Newman, in ; 
Proc. Intern. Conf. Mag., Nottingham (1964). _ - 
» D H. Ridgley and J. H. Geisman, J. Appl. Phya. 39, 592, 

(1 2fS Mo and T. Suzuki, J. Phys. Japan 24, 964 (1968). J 
»R. U Rosenberg, W H Ctoud F. J. Darnell, and R. B, 
Flippen, Phys. Letters 25A, 723 (1967). mi fi968) v 

Roger and R. Fruchart, Mat. Res. Bull. 3, 253 
he. Fwchart, A. M. Triquet, and R. Fruchart, Ann. Chu% 

''.".'pSart, A. M. Blanc, E. Fruchart, J. ^^CNR* 
-J. P. Senateur, CoUoque International No. 157 du C.N.R.* 

(1 " 6 a'm'. Blanc, E. Fruchart, and R. Fruchart, Ann. Chun. | 
251 (1967). 



M 4^T ( Schi and K. Motfeuki, J. Phys. Soc. Japan 24, 742 

(l 2 6 R ) J Gambino, T. P. McGuire, and Y. Nakamura, J. Appl. 

"SSftiS PhTliag. 17, 623 (1968) 
« M. C. Cadeville, these Strasbourg 5 > - P. 29 - 
» A Gerard, Bull. Soc. Belg. Phys. 5, 43 (1966). 
"^Bonnerot, R. Fruchart, and A. Roger, Phys. Letters 26A, 

S3 » RJ^Bailey and J. EDuncan, Inorg. Chem. 6, 1444 (1966) . 
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( -M.C. Cadeville, E. Damel,4nd A J. P .Meyer, Colloq| 
International No. 157 du CJ^.RiS. (1965) , P. 361. 
r- »H. Nowotny and E. Henglein, Monats. Chem.79, 385 (l* 
»iS. Rundqvbt and F. Jellinek, Acta Chem. Scand. 13, 

(1 So!'Antad' and H. Nowotny, Z. Physik. Chem. B38, | 

(1 » 3 H.'Nowotny and E. Henglein, Z. Physik.! Chin... 
(1938). 



BEST AVAILABLE COF^ 



i 

tCH 1969 



sphorous in 
lectrons are 
stitutions of 
ges of mag- 
se found 28 * 29 

!c properties 
type M 2 P. 
brief struo 
; Af 2 P phos- 
r e the same 
d Ni 2 P, are 

.. Newman, in 

Phys. 39, 592 

(1968). 

ell, and R. B. 

J, 253 (1968). 
t, Ann. Chim. 

Bouchaud, and 
dn C.N.R.S. 

Ann. Chim. 2, 

(unpublished) . 
Rend. 255, 339 

eyer, Colloque 

79, 385(1948). 
Scand. 13, 425 

lem. B38, 356 

hem. B40, 281 



t CRYSTALLOGRAPHIC AND MAGNETIC PROPERTIES 



isomorphous with hexagonal Fe 2 P, 31,32 and others have 
the orthorhombic structure 34 of C02P. The structure of 
Cr 2 P has not yet been determined. Rundquist 86 has 
discussed the similarities between the hexagonal and 
orthorhombic phosphides. In order to relate magnetic 
properties to crystal structure, we review briefly this 
similarity here, emphasizing the metalloid environ- 
ment of the metal atoms and the elementary structural 
unit common to both structures. 

Along the [001] direction of the hexagonal structure 
or the [010] direction of the orthorhombic structure, 
phosphorous atoms form, alternately, tetrahedra and 
square-based pyramids, as shown in Fig. 1 (a) , making 
a canal of triangular cross section. Occupation of both 
tetrahedral and pyramidal sites cannot occur simul- 
taneously along the same canal. Consequently, two 
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2. Arrangement of rhombohedral subcells in hexagonal 
FesP structure. 
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Fig. 1. Building blocks of hexagonal and orthorhombic M^P 
structures, (a) Phosphorous canals, (b) Relationship of pyr- 
ramidal-site and tetrahedral-site atoms forming a rhombohedral 
subcell. 

kinds of canals are possible: one has only tetrahedral 
sites occupied, the other only pyramidal sites occupied, 
in both the hexagonal Fe 2 P and orthorhombic C02P 
structures, canals of different type share common 
pyramidal square bases to form canal pairs. This 
pairing creates a sequence of elementary structural units 
along the paired canals. As shown in Fig. 1 (b) , these 
units are rhombohedral and consist of an occupied 
tetrahedron-occupied pyramid pair, and an empty- 
tetrahedron-empty-pyramid pair. The hexagonal Fe 2 P 
structure of Fig. 2 corresponds to an hexagonal assem- 
bly of these canal pairs whereas the orthorhombic Co 2 P 
structure of Fig. 3 results from a zig-zag stacking. 

II. SAMPLE PREPARATION 

Solid solutions were prepared at temperatures be- 
tween 850° and 900°C in evacuated, sealed silica am- 

M H. Nowotny, Z. Anorg. Chem. 254, 31 (1947). 
* S. Rundqvist, Arkiv fur Kemi 20, 7 (1962). 



pules by solid-phase diffusion of phosphides prepared 
by a previously described method. 9 

III. EXPERIMENTAL METHODS 

1. X-ray diffraction: Diffraction patterns were taken 
with a Seeman-Bohlin camera with monochromatic 
chromium Kai radiation. The precision of the meas- 
urements is about 10~* A. 

2. Magnetic analysis: Curie and magnetic- transition 
temperatures were determined with a thermomagnetic 
balance from 77°-1400°K in fields of a few hundred 
Oersteds. Magnetic moments were obtained in fields 
up to 26.6 kOe and in temperatures down to 20.4°K. 

3. Mbssbauer spectra: The apparatus includes a Co 57 
source imbedded in palladium and placed in a parabolic- 
movement vibrator having a parabolic displacement. 
The spectra, obtained from 90° to 650°K, were recorded 
in a four-hundred-channel analyzer. Zone-melted iron 
was used as a standard. 
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Fic. 3. Arrangement of rhombohedral subcells in orthorhombic 
C02P structure. ; 
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6,700 



0.04 for preparation at 600° and 9S0°C, respectively. 
Vegard's law is not obeyed, and d vs x shows an inflec- 
tion point near *=0.60. To obtain a unique phase at 
850°C, it was necessary to add an excess of 3 at.% 
phosphorous throughout the range *<0.50. 

2. (Coi-^Mn») 2 P: Hexagonal structures in this sys- 
tem are restricted to the compositional range *>0.80 
(Fig. 6). Particularly marked anomalies in a! and V 
occur at the composition MnCoP. 

3. (Mni^Cr*) 2 P: Substitution of 10 at.% chromium 
into Mn 2 P is sufficient to make the hexagonal struc- 
ture disappear. For s>0.10, the system is multiphase. 

4. (Coi_JSTix)2P: In this system, hexagonal struc- 
tures are found over most of the compositional range, 
orthorhombic structures being restricted to the range 
s<0.15. The hexagonal c parameter has a broad mini- 
mum at #=0.65. 

C. Complete Solid Solutions 

The hexagonal system (Fe^NWjP has complete 
solid solutions. The lattice parameters vary anomalously 
with x (Fig. 7), the most accentuated changes occurring 
at the composition x=0.8S, which has the minimum 
volume. 



Fig 4 Lattice-parameters and structurallhomogeneity ranges 
for the system (Fei^MnJaP. 

IV. CRYSTALLOGRAPHIC STUDY 
A. Systems with Three Single-Phase Regions 

1. (Fei-xMn x ) 2 P: The solid-solution system (Fe t _^- 
Mn^P has orthorhombic symmetry in the inter- 
mediate compositional range 0.31<s<0.62 (Fig. 4). 
The orthorhombic and hexagonal phases are separated 
by very narrow two-phase regions, and the limiting 
compositions are very slightly dependent on the tem- 
perature of sample preparation. Lattice parameters do 
not obey Vegard's law, and abrupt changes in a! occur 
aUs= 0 >5} or MnFeP. The hexagonal phase is more dense 
than the orthorhombic phase (107.58 A 3 vs 108.19 A 

for*=0.62.) m m 

2. (Mni_JSfix)aP: This system is similarly character- 
ized by orthrhombic structures in the intermediate 
compositional range 0.20<x<0.50. 

3. (Fe^Cr x ) 2 P: This system has orthorhombic 
structures in the compositional range 0.15<s<0.55, 
and the rate of change with x of all three orthrhombic 
parameters a', c\ show anomalous changes at *= 
0.50. However, for s>0.60 the system is multiphase. 

B. Systems with Two Single-Phase Regions 

1 (Coi-JFe,) 2 P: In this system hexagonal struc- 
tures are restricted to the relatively narrow composi- 
tional range s>0.84 (Fig. 5). Here the limiting com- 
positions of the very narrow two-phase region, of 
width A*, vary with thermal treatment: Aa«0.01, 



V. MAGNETIC STUDY 

1. (Coi-JVln^P: Although the binary compounds 
Mn 2 P and Co 2 P are not ferromagnetic, the orthorhombic 
domain of the system (Co!_*Mn x ) 2 P contains a large 
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Fig. 5. Lattice parameters^and structural homogeneity 
for the system (Coi_«Fe»)iP. . 
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Fig. 6. Latticejparameters^and structural homogeneity ranges 
for the system (Coi_,Mn x ) 2 P. 

ferromagnetic compositional region (Fig. 8). The 
Curie temperatures T c vary rapidly with x, exhibiting a 
sharp maximum of 583°K at a; = 0.5 (CoMnP). In the 
interval 0.56 <#< 0.66, there is a metamagnetic to 
ferromagnetic transition with increasing temperature. 
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Fig. 8. Magnetic transition'temperatures'and magnetization at 
0°K for the system (Coi-*Mn») t P. 

The transition temperature T\ increases rapidly with 
x while T c decreases, the domain of ferromagnetism in 
weak external fields terminating at s=0.66. T\ varies 
strongly with applied field (Figs. 9-11) , and for x>0.66 
the magnetic order-disorder transition is a N6el tem- 
perature T N (Figs. 8-10). The saturation magnetiza- 
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Fig. 9. Thermomagnetic curves at various applied field strengths 
(Co 0 . w Mn fc#4 ) t P.7^ 
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Fig. 10. Thermomagnetic curves at various applied field strengths 
for (Coo.i4Mno.«)iP. 
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Fig 11. Magnetization curves at 20.4°K for various compositions 
(Coi-xMn^iP that are metamagnetic at low temperatures. 
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Fig. 13. Curie temperature for the system (Fei_*Ni*)iP. 

tion, extrapolated to 0°K, is also shown in Fig. 8. It 
reaches a maximum of 3.06 /i B /molecule at s«0.55. 

2. (Coi_xFe x )2P: As in the previous system, the 
Curie temperature T c rises rapidly with x to a maximum 
of 459°K at s=0.62 (Fig. 12). However, here there is 
no metamagnetic state, and there is a discontinuity in 
T c at x=0.84, where the structural change occurs, from 
407°-441°K. Since thermomagnetic analysis of samples 
with x<0.50 revealed the presence of a magnetic im- 
purity (the amount increasing with preparation tem- 
perature) having a Curie temperature slightly less 
than that of pure cobalt, for s<0.50, the phosphorous 
content was systematically increased by 3 at.% in these 
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Fic 2 12. Curie temperature and magnetization at 0«K for the Fig. 14. Magnetization vs applied field at 20.4'K of Feigg 



system (CowFe«)iP. 
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compositions in order to obtain single magnetic phases. 
The magnetization (Fig. 12) decreases abruptly from 
2.75 to 1.83 /x B /molecule on passing from the hexagonal 
to the orthorhombic phase, and within the orthorhom- 
bic phase it reaches a maximum of 2.02 /is/molecule 
at s=0.60. 

3. (Fe^Nix^P: Substitution of Ni into Fe 2 P pro- 
duces a similar sharp increase in T e (Fig. 13), which 
reaches a maximum of 342°K at *=0.08, decreasing 
rapidly with larger x to less than 100°K at z=0.50. 
The saturation magnetization has not yet been meas- 
ured. 

4. (Fei_ x Mn 3C ) 2 P: Very small substitutions of man- 
ganese into Fe 2 P induce metamagnetism. By #=0.03, 
the antiferromagnetic state is stable in fields up to 
26.6 kOe (see Fig. 14). Similar changes occur with 
chromium substitutions. 

5. Fe 2 _eP: Cadeville and Meyer 15 have reported 
variations in magnetization with deviations from 
stoichiometry as large as 2.76 At B /molecule for Fe 2 P to 
2.14 MB/molecule for Fe 2 _«P. Similarly, we find that 
two single-phase samples of nominal compositions 
(Fe 2 P), but differing by one weight percent iron, have 
magnetizations of 2.67 /z B and 2.44 MB/molecule. 

VL MOSSBAUER SPECTRA 
A. Studies in the Paramagnetic State 

The Mossbauer spectrum of Fe 2 P has already been 
studied. 16 - 18 Pyramidal-site iron atoms are subject to a 
significant electric-field gradient that splits the reso- 
nance line by ^0.4 mm/sec. The tetrahedral-site iron 
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Fig. 16. Mossbauer spectra of (Fe,Coi_*)iP in the 
paramagnetic state. 
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Fig. 15. Mossbauer spectra of (FcCri-JjP in the 



atoms exhibit a single resonance line of enlarged width, 
probably reflecting a slight deformation of the tetra- 
hedron. Therefore, Mossbauer studies of (Fei_xAf x ) 2 P 
systems permit a direct determination of the site pref- 
erences of the solute M atoms. 

1. M =Mn or Cr: Changes with x in the Mossbauer 
spectra indicate a total substitution of manganese and 
chromium atoms for pyramidal-site iron. Only the 
single line with enlarged width remains in the spectra 
of FeMnP and FeCrP (Fig. 15) . 

2. M = Co or Ni: Cobalt and nickel substitute prefer- 
entially for tetrahedral-site iron. About 15% of the 
cobalt atoms occupy pyramidal sites in FeCoP. For 
z>0.70, only the symmetric doublet, separated by 
0.83 mm/sec, appears in the spectrum (Fig. 16). About 
25% of the nickel atoms occupy pyramidal sites in 
FeNiP. 

B. Studies in the Ferromagnetic State 

The spectra of the system (Fei^ c Co I ) 2 P at T= 
294°K<r e for a series of compositions in the range 
0<z<0.25 vary regularly with x and show little change 
on passing from hexagonal to orthorhombic symmetry. 
They have been interpreted with the method of 
Wertheim, 36 where the parameters to be determined are 
the internal fields at the tetrahedral-site and pyramidal- 
site nuclei F„' and H n p , the quadrupolar field strengths 
ieqV aa at pyramidal sites (those at tetrahedral sites 
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Table I. Mossbauer parameters at 294°K for 
(Coi^Fe,)»P- 



X 


(kOe) 


(kOe) 




1 


e 


0.8 


67 


132 


0.63 


0.46 


90° 


0.75 


69 


141 


0.67 


0.56 


90° 


0.70 


69 


153 


0.67 




20° 


0.50 


69 


153 


0.74 




20° 



assumed negligible) , an asymmetry parameter 17, and 
the angle 8 between the spin direction and the principal 
axis of the electric-field gradient in the pyramidal sites. 
The results are summarized in Table I. By symmetry, 
the principal axis of the electric-field gradient should 
be perpendicular to the phosphorous canals, so that 
a 0=90° means that the spins are parallel to the c 
axis in the hexagonal structure, to the V axis in the 
orthorhombic structure. 

In the range s>0.25, an increase in the complexity of 
the spectra (Fig. 17) indicates a spin reorientation. In 
first approximation and neglecting the asymmetry 
parameters, a 0^20° is obtained with the aid of Kun- 
dig's 37 graphs. Since the axes of the pyramidal sites make 
an angle of about 30° with the orthorhombic c' axis, 
it is reasonable to assume that the spins have become 
oriented parallel to this axis. For #=0.25, the spins 
rotate continuously with temperature in the V-c' 




V {mm /uJ 
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Fig, 17. Rotation of spins as a function of composition 
in (Fe.Co,^)^.' 



plane: at 90°K they are parallel to c' and at 294°K to 
b' (Fig. 18). These results are confirmed by x-ray 
diffraction on powder samples oriented in a magnetic 
field at 294°K. A detailed, complete analysis of the 
spectra will be published soon. 38 

vn. DISCUSSION 

Mixed phosphides of formula M 2 P have a remarkable 
alternation of hexagonal and orthorhombic structures 
with increasing electron/atom ratio, 

Mn 2 P(hex) -MnFeP (ortho) -Fe 2 P (hex) 

-CosP (ortho) -Ni 2 P(hex) . 

Anomalies in the variation with composition of 
lattice parameters and magnetic properties, which are" 




Fig. 18. Rotation of spins as a function of temperature 
in (Fe«. n Coo.«)iP. 



particularly pronounced in the system (Coi-JMn*)sP, 
can be attributed, as a result of the Mossbauer studies, 
to the site preferences of the solute atoms. The com- 
pletely ordered compounds FeMnP, FeCrP, MnCoP, 
and, without doubt, MnNiP, belong to the important 
family of orthorhombic E phases represented by 
ZrFeP N' 40 ; a homologous family of ordered hexagonal 
phases would be represented by NbMnSi. 41 ' 42 For many 
members of these families, the ordering energy may 
be so large that the homogeneity range is,strictly limited 
to the ideal compositions, as in ZrFeP and NbMnSi. 
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Not only the substitutional ordering, but also deviations 
from stoichiometry of the metal and not the metalloid 
atom, 41 suggest that the two structures Fe 2 P and Co 2 P 
should both be viewed as metalloid skeletal arrays 
having twice as many interstitial metal atoms as 
metalloid atoms and the common, rhombohedral subcell 
of Fig. 1(b). 

Although incomplete, magnetic studies reveal a 
fundamental difference between the phosphides Af 2 P 
and MiV: In the systems (Fei_ x Mx) 3 P, where M = Co 
or Ni, the linear change with x in magnetization 13 ' 26 - 27 
agrees with the hypothesis of filling of incomplete 3d 
half-bands, cobalt contributing one and nickel two 
extra electrons. The magnetization changes by 0.94 
p B /cobalt atom and by 1.92 ^B/nickel atom. 27 Likewise, 
the abrupt decrease with x of the magnetization if 
M=Cv shows that the scheme proposed 29 for the 
borides M 2 B and MB applies equally to the M 3 P 
phosphides. This scheme uses the virtual-state concept 
of Friedel, 44 ' 45 which was introduced to account for the 
magnetic properties of the alloys of the first-row tran- 
sition metals. Thus the behavior of the compounds 
like M 2 B and MB, is essentially the same as for 
the constituent metals, except for the contribution of 
the metalloid electrons to the 3d bands. Here the metal- 
atom array forms the skeletal net. On the other hand, 
the phosphides Af 2 P, although poor in metalloid atoms, 
exhibit a behavior reminescent of boron- or phosphorous- 
rich compounds and compounds having a pronounced 

«S. Rundqvist, Acta Chem. Scand. 14, 1961 (1960). 
44 J. Friedel, Nuevo Cimento Sup. VII 2, 287 (1958). 
« J. Friedel, J. Phys. Radium 23, 692 (1962). 



ionic character, such as oxides and sulfides, for which 
the anionic array provides the structural skeleton. This 
curious contrast is supported by two additional obser- 
vations: 

1. Although maximum Curie temperatures are high 
(583°K), nevertheless ferromagnetism is confined to a 
relatively narrow range of electron/atom ratios in 
both the (Coi_*Mn x ) 2 P system, which is completely 
ordered, and the (Coi_ <c Fe z ) 2 P system, which is partially 
ordered. (The maximum T € would probably occur for 
CoFeP were the atomic ordering complete.) 

2. Ferromagnetism is replaced by metamagnetism 
for electron/atom ratios a little smaller than that of 
CoMnP in the system (Coj_JV!n x ) 2 P. Therefore, it is 
reasonable to anticipate that Fe 2 P, which should have a 
similar 3d electron/atom ratio, exhibits a magnetiza- 
tion that is extremely sensitive to any reduction in this 
ratio (whether by substitution of Mn or Cr or by devia- 
tion from stoichiometry, Fe^P) because it is at the 
critical ratio at which ferromagnetism gives place to 
metamagnetism. Although not well established by 
experiment, it appears from the Mossbauer spectra, 38 
that small reductions in the electron/atom ratio induce 
progressively a partial canting of the spins in Fe^. 
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A series of iron, cobalt, and nickel metal phosphides of chemical 
formula Fe 2 P, CoP, and Ni 2 P with specific surface areas of around 
3 m 2 g" 1 were synthesized by means of temperature-programmed 
reduction (TPR) of the corresponding phosphates. These phos- 
phides were also successfully prepared in dispersed form on a 
silica support (90 m 2 g -1 ) for use as catalysts. The phase purity 
of these materials was established by X-ray diffraction (XRD), 
and surface properties were determined by N 2 BET specific sur- 
face area (5 g ) measurements and CO uptake determinations. The 
activity of the silica-supported catalysts in hydrodenitrogenation 
(HDN) and hydrodesulfurization (HDS) was evaluated in a three- 
phase trickle-bed reactor using a model liquid feed containing 
2000 ppm nitrogen as quinoline, 3000 ppm sulfur as dibenzoth- 
iophene, 500 ppm oxygen as benzofuran, 20 wt% aromatics as 
tetralin, and balance aliphatics as tetradecane. The reactivity study 
showed that the HDS activity sequence for the three samples 
was Ni 2 P/Si0 2 > CoP/Si0 2 > Fe 2 P/Si0 2 , while the HDN activity 
followed the sequence CoP/Si0 2 > Ni 2 P/Si0 2 > Fe 2 P/Si0 2 . Com- 
pared with a commercial Ni-Mo-S/y-Al 2 0 3 catalyst, Ni 2 P/Si0 2 
had a higher HDS activity (90 vs 76%), but a lower HDN activ- 
ity (14 vs 38%), based on equal sites loaded in the reactor. The 
sites were determined by CO chemisorption for the phosphide 
and low-temperature 0 2 chemisorption for the sulfide. XRD and 
X-ray photoelectron spectroscopy characterizations of the spent 
catalysts indicated that the Ni 2 P/Si0 2 catalyst was tolerant of 

SUlfur. © 2002 Elsevier Science (USA) 

Key Words: hydrodesulfurization; hydrodenitrogenation; transi- 
tion metal phosphides; Fe 2 P; CoP; Ni 2 P. 



INTRODUCTION 

Demands for a cleaner environment have led to a global 
tightening in the allowed sulfur content in fuels and in- 
creased restrictions on the release of nitrogen oxides. For 
example, in the case of sulfur the U.S. Environmental Pro- 
tection Agency (EPA) has issued regulations that would 
lower its allowed content in diesel fuel from the current 500 

1 Present address: Luna Innovations, 2851 Commerce St., Blacksburg, 
VA 24073. 

2 To whom correspondence should be addressed. E-mail: oyama@ 
vt.edu. 




to 15 ppmw in 2006, and in gasoline from 300 to 30 ppmw 
by 2004 (1, 2). For this reason there are considerable ef- 
forts being expended to develop new technologies for the 
production of clean fuels, like adsorption, extraction, oxida- 
tion, alkylation, and bioprocessing (3). Currently, however, 
hydroprocessing appears to be the technologically prefered 
solution (3). Hydroprocessing refers to a variety of catalytic 
hydrogenation processes that saturate heteroatomatic rings 
and remove S, N, O, and metals from different petroleum 
streams in a refinery (4). Because of the tighter environ- 
mental regulations new types of catalysts, which are eco- 
nomic, have long life, and possess high activity, are highly 
desired. In this work we present results on a new type of 
hydroprocessing catalyst: transition metal phosphides. 

Transition metal phosphides have attracted considerable 
interest for some time because these materials are tech- 
nologically important as semiconductors, luminescent de- 
vices, and electronic components (5). A brief thermody- 
namic analysis of the potential stability in H 2 S was carried 
out early in 1975 (6), and it was revealed that this group 
of materials is potentially stable and sulfur resistant. How- 
ever, the transition metal phosphides as a class of materials 
have received little attention in the field of catalysis. This 
is probably because the development of synthetic meth- 
ods of producing materials with high surface areas was not 
achieved. 

The combination of the iron group metals Co and Ni 
with Mo and W in commercial hydroprocessing catalysts 
(7-9) and the use of phosphorus as a promoter (10-12) 
is well-known. Many workers have studied the effect of 
phosphorus in sulfide catalysts (13-16), and the topic has 
been reviewed by Iwamoto and Grimblot (17). It was con- 
cluded that in these materials, the phosphorus was found as 
a phosphate and primarily modified the properties of the 
I support, and only indirectly the active phase. For example, 
I phosphorus altered the acid-base character of alumina and 
improved dispersion of molybdenum on the support. It also 
enhanced the solubility of the precursor metals in the prepa- 
' ration stages and allowed the synthesis of high-loading cata- 
> lysts (18). The effect of phosphorus strongly depended on 
its content, with an effect that was usually negative at high 
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loadings. Phosphorus showed no effect or a small positive 
effect on the hydrodesulfurization (HDS) of thiophene, and 
a positive effect on the hydrodenitrogenation (HDN) of 
quinoline (19), pyridine (20), and piperidine (11). Notwith- 
standing the considerable studies on promoter effects, until 
recently phosphorus compounds in the form of phosphides 
had not been examined in the hydroprocessing field. The 
first report of the use of the iron group metal phosphides as 
hydrodenitrogenation catalysts was by Robinson et al (21), 
who prepared Co 2 P and Ni 2 P on silica, alumina, and car- 
bon. They reported that carbon and silica were the best sup- 
ports and that Ni 2 P, in particular, was very effective in HDN. 
These studies were carried out at high conversion, and an 
assessment of the intrinsic activity is not possible. The first 
studies on Mo and W phosphides were carried out by Li 
et al (22), Oyama et al (23), and Clark et al (24), who re- 
ported good activity for HDS and HDN. Studies in the Prins 
group (25) confirmed the activity in HDN, and in a recent 
comparison of Co 2 P, Ni 2 P, MoP, CoMoP, and NiMoP (26) it 
was concluded that the areal activity of MoP was the high- 
est. Earlier, the olefin hydrogenation activity of Ni 2 P sup- 
ported on alumina and other phosphides was explored by 
Nozaki and coworkers (27-29). It was found that the hydro- 
genation activity for butadiene drastically decreased, in the 
order Ni 2 P > Co 2 P > FeP. They also reported that a trace of 
oxygen could increase the activity of Ni 2 P while lowering 
the activity of Ni for the butadiene hydrogenation reaction. 
Nickel-phosphorus alloys have been reported also in amor- 
phous form and their activities for hydrogenation have also 
been studied (30, 31). The amorphous alloy was prepared 
by an electroless plating technique from mixtures of sodium 
citrate, nickel sulfate, sodium hypophosphite, sodium ac- 
etate, and a silica gel support, or by the chemical reduction 
of nickel acetate and sodium phosphate with sodium boro- 
hydride. The supported material was subjected to various 
treatments, including oxidation at 403 K and reduction in 
H 2 at 553 K. The catalyst was found to be active for the hy- 
drogenation of nitrobenzene (32) and benzaldehyde (33), 
with a turnover rate similar to that of Ni (1 x 10" 3 s" 1 ). 
In summary, although there are some reports concerning 
the catalytic behavior of these phosphides, few studies have 
concentrated on the subject of hydrodesulfurization and hy- 
drodenitrogenation for application in petroleum refining. 

The present work presents an in-depth study of the 
preparation of transition metal phosphides of the iron 
group (Fe, Co, Ni) and their evaluation in the hydropro- 
cessing of a model feed mixture. Initially unsupported 
bulk materials were prepared to provide a reference for 
the synthesis of supported materials. Several metal-to- 
phosphorus (M/P) ratios were explored to ascertain which 
stable phases could be prepared and to determine the 
conditions for temperature-programmed reduction. Sub- 
sequently, the preparation was extended to the supported 
system to obtain materials of high surface area suitable for 
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catalytic testing. Silica was chosen as the carrier to mini- 
mize support effects and make possible the elucidation of 
the intrinsic catalytic activity of the phosphides. 

EXPERIMENTAL 

Materials 

The support used in this study was a fumed silica (Cabosil, 
L90). The precursors for iron, cobalt, and nickel were 
Fe(N0 3 ) 3 • 9H 2 0 (Aldrich, 99.99%), Co(N0 3 ) 2 ■ 6H 2 0 
(Aldrich, 99.99%), and Ni(N0 3 ) 2 • 6H z O (Aesar, 99%), re- 
spectively, while the precursor for P was ammonium or- 
thophosphate (NH4) 2 HP0 4 (Aldrich, 99%). The chemi- 
cals utilized in the reactivity study were dibenzothiophene 
(Aldrich, 99.5%), quinoline (Aldrich, 99.9%), benzofuran 
(Aldrich, 99.9%), tetralin (Aldrich, 99.5%), and tetrade- 
cane (Jansen Chimica, 99%). The gases employed were He 
(Airco, Grade 5), CO (Linde Research Grade, 99.97%), 
0.5% 0 2 /He (Airco, UHP Grade), H 2 (Airco, Grade 5), N 2 
(Airco, 99.99%), and 30% N 2 /He (Airco, UHP Grade). 

Synthesis 

Unsupported bulk transition metal phosphides were pre- 
pared in two steps. In the first step, phosphate precursors 
were synthesized by reacting metal nitrates with ammo- 
nium phosphate, and in the second step, these phosphates 
were reduced to phosphides by the method of temperature- 
programmed reduction. Because the procedures for prepar- 
ing unsupported bulk iron, cobalt, and nickel phosphates 
were similar, the preparation of iron phosphate (FeP0 4 ) is 
used here to illustrate the process. First, 4.9 g (37.13 mmol) 
of ammonium phosphate (NH 4 ) 2 HP0 4 was dissolved in 
300 cm 3 of distilled water to form a transparent colorless 
solution, and 15 g (37.13 mmol) of iron nitrate (Fe(N0 3 ) 3 • 
9H 2 0) was then added. The clear solution immediately 
turned into a light color mixture with some precipitate, 
but stirring resulted in the formation of a transparent so- 
lution. In the case of nickel and cobalt, several drops of 
nitric acid were needed to give rise to a homogenous so- 
lution. The water was then vaporized from the solution 
on a hot plate and the resulting paste was dried at 393 K 
for 3 h and calcined at 773 K for 6 h in an oven. The 
amount collected (with some minor losses) was 5.54 g, 
which corresponds to 36 mmol of iron phosphate of formula 
FeP0 4 • H 2 0. The phosphate was then ground with a mor- 
tar and pestle and sieved to 16/20 mesh (0.65- to 1.2-mm- 
diameter particles). In the second step of preparation, 
temperature-programmed reduction (TPR) was utilized 
to convert the phosphate into phosphide. The reduction was 
carried out in a U-shaped quartz reactor placed in a furnace 
controlled by a temperature programmer (Omega Model 
CN 2000). The temperature was raised at p = 0.0167 K s" 1 
(1 Kmin -1 ) and was monitored by a local chromel-alumel 
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TABLE 1 

Quantities Used in the Preparation of Supported Samples 





Materials used 


Catalyst properties 




Metal 




Phosphide Metal 




nitrate 


(NH 4 ) 2 HP0 4 


loading loading 


Sample 


Silica (g) (mol) 


(mol) 


(wt%M,P) (mol%M) 


Fe 2 P/Si0 2 


20 0.0462 


0.0231 


14 11 


C0P/S1O2 


20 0.0231 


0.0231 


9.4 6.2 


Ni 2 P/Si0 2 


20 0.0231 


0.0231 


9.4 6.1 



thermocouple placed in a thermowell near the center of 
the reactor bed. The H 2 flow rate was set at 1000 /-tmol s" 1 
(1500 cm 3 min" 1 ) per g of sample. A portion of the exit gas 
flow was sampled through a leak valve into a mass spec- 
trometer (Ametek/Dycor Model MA 100) and the masses 
2(H 2 ), 4(He), 18(H 2 0), 28(N 2 ), 32(0 2 ), 34(PH 3 ), 15(NH), 
44(C0 2 ), 31(P), and 62(P 2 ) were monitored during the 
experiment, and these were recorded together with the tem- 
perature by an online computer. At the end of the temp- 
erature program, the sample was cooled in helium to room 
temperature and was passivated in a 0.5% 0 2 /He flow 
for 2 h. 

Phosphide samples supported on a silica support were 
prepared by modifying the two-step procedure used in the 
synthesis of the unsupported samples. Aqueous phosphate 
solutions were obtained as before and were used to impreg- 
nate silica by the incipient wetness impregnation method. 
The quantities used for the three supported samples are 
listed in Table 1. Prior to use, the silica was dried at 393 K 
for 3 h and calcined at 773 K for 6 h and was found to have 
an incipient wetness point of 2.2 cm 3 g" 1 . After impregna- 
tion, the powders were dried at 393 K for 3 h and calcined at 
773 K for 6 h. The calcined samples were ground with a mor- 
tar and pestle, pelietized with a press (Carver, Model C), 
and sieved to 16/20 mesh size. The TPR process was similar 
to that used for the bulk samples using the same heating 
rate, p = 0.0167 Ks" 1 (IK min" 1 ). 

Characterization 

The synthesized materials were characterized by CO 
chemisorption, N 2 physisorption, and X-ray diffraction 
(XRD) measurements. Irreversible CO uptake measure- 
ments were used to titrate the surface metal atoms and to 
provide an estimate of the number of active sites on the 
catalysts. Uptakes were obtained after passivated samples 
were rereduced and are denoted ex situ in this paper. Usu- 
ally, 0.2 g of sample was loaded into a quartz reactor and 
treated in H 2 at 723 K for 2 h. After cooling to room tem- 
perature in He, pulses of CO in a He carrier flowing at 
27 /xmol s" 1 (40 cm 3 NTP min -1 ) were injected through a 
sampling valve and the 28 (CO) signal was monitored with 
a mass spectrometer. Uptakes of CO were also measured 



for the spent samples. The procedures were the same as 
those used for the fresh samples. Prior to the measurement, 
the spent samples removed from the hydrotreating reactors 
were washed in hexane and dried. 

BET surface area measurements were carried out right 
after the CO uptake determinations, using a similar tech- 
nique. Adsorption at liquid nitrogen temperature was per- 
formed using a 30% N 2 /He stream, and the desorption area 
obtained after rapid heating was compared to the area of 
a calibrated volume (35.4 /imol). The surface area was cal- 
culated from the one-point BET equation, which is reason- 
able for nonmicroporous materials such as those used here. 
X-ray diffraction (XRD) patterns of the samples were de- 
termined with a Scintag XDS-2000 powder diffractometer 
operated at 45 kV and 40 mA, using Cu Ka monochro- 
matized radiation (A. = 0.154178 nm). The crystallite size 
of the supported sample was calculated using the Scherrer 
equation, D c = Kk/pcos(0), where K is a constant taken 
as 0.9, k is the wavelength of the X-ray radiation, p is the 
width of the peak at half-maximum, corrected for instru- 
mental broadening (0.1°), and 20 is the Bragg angle (34- 
36). The near-surface composition of the nickel samples 
was obtained by X-ray photoelectron spectroscopy (XPS) 
(Perkin-Elmer, Model 5300 with a Mg source) operated 
at 15 kV and 30 mA. The 285.0-eV binding energy peak 
of adventious carbon was used as reference. In the case of 
the spent catalysts, samples were removed from the reac- 
tor, washed in hexane, heated in H 2 to 673 K, and then 
passivated. Since the samples were exposed to the atmo- 
sphere and not sputtered, contamination by carbon from 
the atmosphere was present. 

Reactivity Studies 

Hydrotreating activities of the samples were obtained 
in a three-phase trickle-bed reactor for hydrodenitrogena- 
tion (HDN) and hydrodesulfurization (HDS) with a model 
petroleum liquid containing 2000 ppm nitrogen (quinoline), 
3000 ppm sulfur (dibenzothiophene), 500 ppm oxygen (ben- 
zofuran), 20 wt% aromatics (tetralin), and balance aliphat- 
ics (tetradecane). The operating conditions were close to 
industrial conditions, 3.1 MPa and 643 K, with a liquid 
flow rate of 5 cm 3 /h and a hydrogen flow rate of 100 /xmol 
s _1 (150 cm 3 min* 1 ) corresponding to a gas-liquid ratio of 
9800 SCF H 2 /barrel. The detailed description of the testing 
system is reported elsewhere (37). Quantities of catalysts 
loaded in the reactor correspond to the same amount of 
ex situ CO uptake (35 /zmol). Prior to reactivity measure- 
ments, the catalyst samples were pretreated in exactly the 
same manner as before the ex situ CO uptake determina- 
tions. Hydrotreating samples were collected every 2 or 3 h 
in sealed septum vials and were analyzed off-line with a gas 
chromatograph (Hewlett Packard, 5890 A) equipped with a 
0.32 mm i.d. x 50 m fused silica capillary column (CPSIL- 
5CB, Chrompack, Inc.) and a flame ionization detector. 
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RESULTS AND DISCUSSION 

Properties and Preparation of Bulk Phosphides 

The crystal structures of the Fe, Co, and Ni phosphides are 
shown in Fig. 1 and their lattice parameters are summarized 
in Table 2. Both Fe 2 P (38) and Ni 2 P (39) adopt the same 
hexagonal structure (space group: Pfrm), while CoP (40) 
takes on an orthorhombic structure (space group: P nm a). 
The crystal structures (Fig. 1) and lattice parameters can be 
used to calculate the bulk density (p) and the surface metal 
atom density (h) of the solids (Table 2). For the Fe 2 P and 
Ni 2 P samples there are two, three, and two atoms on the ac, 
ab, and be unit cell faces, respectively. For CoP, every unit 
cell face has two atoms. 

The investigation of the phosphide materials in this study 
was begun by a study of the synthesis of bulk materi- 
als of various metal-to-phosphorus ratios (M/P) by the 
temperature-programmed method. This was carried out to 
prepare suitable references for the supported materials, and 
to evaluate the conditions of preparation needed when a 
support was employed. 

The synthesis of the bulk phosphides involved two stages, 
preparation of phosphate precursors and reduction of the 
precursors in a temperature-programmed manner. The re- 
sults for iron are discussed first. The precursors prepared 
in the first stage had different Fe/P ratios, set at the time of 
preparation by adjusting the mole ratio of the constituents. 
The precursors were formed by the thermal decomposition 
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FIG. 2. Temperature-programmed reduction of bulk iron phosphates 
at p = 1 K/min (0.01667 K s" 1 ) (Fe-to-P ratios used in the preparation are 
indicated). 



of mixtures of the metal nitrate and ammonium phosphate 
in air. Because the nitrate and ammonium ions are unsta- 
ble at high temperature, it was expected that the mixture 
would decompose to a metal phosphate or metal oxide- 
phosphate. At high metal ratios some metal oxide was also 
likely. X-ray diffraction patterns of the precursors showed 
that they were amorphous. 

Temperature-programmed reduction (TPR) of the pre- 
cursors was carried out in the second stage (Fig. 2). Only 
the results for mass 18(H 2 0) are shown, as the other mon- 
itored masses were featureless or provided little additional 
information. The TPR traces show two systematic trends 
for the samples with different Fe/P ratios. First, for higher 
Fe contents (Figs. 2a-2c), the TPR traces show more peaks 
and a more complicated overall reduction pattern. Second, 
for higher Fe contents, all the reduction peaks including the 
initial and final features are shifted to lower temperatures. 
These trends are reasonable. For higher Fe contents, the 
precursor is probably a mixture of iron oxide, iron phos- 
phate, and possibly other components, and their separate 
reduction results in different peaks. Also for these higher Fe 
contents, the proportion of iron oxide should increase, and 
since this oxide is easily reducible, its TPR peaks appear at 
lower temperatures. In fact, it is likely that some metallic 



TABLE 2 

Surface Metal Density and Bulk Density of Fe 2 P, CoP, and Ni 2 P 



Sample 



Lattice parameter (nm) 



Surface metal density (10 15 atoms cm" 2 ) 
ab plane be plane ac plane Average (n) 



Compound density (p) 
(g cm" 3 ) 



Fe 2 P 
CoP 
Ni 2 P 



0.5867 0.5867 0.3458 1.01 0.986 0.986 0.994 

0.5077 0.3281 0.5587 1.20 1.09 0.706 0.999 

0.5859 0.5859 0.3382 1.01 1.01 1.01 1.01 



6.67 
6.20 
7.09 
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Fe is formed and that it assists in the reduction of the other 
components. For low Fe contents (Fig. 2d), there is no iron 
oxide and reaction occurs at the intrinsic reduction temper- 
ature of the main phosphate phase (—920 K). Because this 
temperature is high, even if a mixture of components other 
than phosphates existed, individual reduction steps cannot 
be resolved and the whole process appears to occur in one 
stage. 

Analysis of the products of TPR was carried out by XRD 
(Fig. 3). The diffraction patterns all show a high background 
because of fluorescence by the iron. However, the pres- 
ence of distinct iron phosphides can be seen easily. For 
the samples with Fe/P ratios of 3/1 and 2/1, the XRD re- 
sults show the expected phases of Fe 3 P (Fig. 2a) and Fe 2 P 
(Fig. 3b). Comparison is made with standards from the pow- 
der diffraction file (PDF) (41), as indicated in the figure. For 
the sample with Fe/P ratio of 1/1, the obtained phase was 
still Fe 2 P (Fig. 3c), and there was a deficiency in phospho- 
rus. Similarly, for the sample with an Fe/P ratio of 1/2, the 
observed phase was FeP. Thus, although the preparations 
were carried out with stoichiometric quantities of metal and 
phosphorus, the final products tended to be metal-rich. 

There are probably several processes which contributed 
to the loss of phosphorus in these samples. Some of the loss 
may have occurred during the TPR process, as traces of PH3 
were detected in the mass spectrometer signal and some 
volatile products were observed to condense at the exit of 
the reactor. Some of the loss probably also occurred in the 
calcination step at 773 K to form the phosphate precursor. 
A small amount of a white solid was found to have sublimed 
onto the lid of the ceramic calcination vessel. This was likely 
to be P2O5, with a melting point of 563 K and a sublimation 
temperature of 787 K. The loss of P was observed in all 
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FIG. 3. X-ray diffraction patterns of iron phosphides (PDF file refer- 
ences are also included). 
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FIG. 4. Temperature-programmed reduction of bulk cobalt phos- 
phates at 0 = 1 K/min (0.01667 K s" 1 ) (Co-to-P ratios used in the prepa- 
ration are indicated). 



samples with high P levels. There is also a possibility that 
some extra phosphorus in amorphous form remained mixed 
in the samples, and that the observed phases were simply 
the stable ones under the preparation conditions. 

The trends in the TPR results for the cobalt samples 
(Fig. 4) are similar to those of the iron samples. The reaction 
traces consist of two main features, a low-temperature peak 
and a more complicated high- temperature signal. (The sam- 
ple with Co/P =1/1 has a small feature at —350 K, probably 
due to dehydration). The low-temperature peak appears 
between 500 and 600 K and is more intense for higher Co 
contents (Figs. 4a-4c). It also shifts to lower temperature 
with increasing Co content. This behavior is consistent with 
the reduction of a cobalt oxide species. The higher temper- 
ature signal appears between 700 and 900 K and probably 
corresponds to the reduction of cobalt phosphate, which is 
expected to be more difficult to reduce. It is the main fea- 
ture for the Co/P = 1/1 sample (Fig. 4d). The XRD patterns 
of the reduced cobalt samples (Fig. 5) show that C02P was 
obtained from the samples with Co/P ratios of 3/1, 2/1, and 
3/2 (Figs. 5a-5c), and CoP was obtained as expected from 
the sample with a Co/P ratio of 1/1 (Fig. 5d). It appears that 
C02P is a particularly stable phase at these conditions. 

The same trends in the TPR traces were also observed 
with the nickel samples (Fig. 6). The reduction features 
were a composite of different peaks, occurring at a low- 
temperature range of 500-700 K, and at a high-temperature 
range of 700-850 K. These various features are attributed 
again to the reduction of different compounds, likely nickel 
oxides and nickel oxide phosphate at the low tempera- 
tures and nickel phosphate at the high temperatures. As 
expected, Ni 3 P was observed in the sample with a Ni/P ra- 
tio of 3/1 (Fig. 7a). But for the Ni/P ratios of 2/1, 3/2, and 1/1 
(Figs. 7b-7d), the phase obtained was Ni 2 R Probably, this is 
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FIG. 5. X-ray diffraction patterns of reduced cobalt phosphates (PDF 
file references are also included). 
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the stable phase at the conditions of the preparation, since 
not all the phosphorus in the precursor phosphate mixtures 
is likely to have sublimed. Different from the other samples, 
the TPR trace of the sample with the Ni/P ratio of 1/1 did 
show some P and PH 3 during the reduction process. This 
accounts for the loss of P in this sample. 

To summarize these TPR and XRD results, the reduction 
of phosphate precursors occurs readily in the Fe, Co, and 
Ni systems, with maximum temperatures of about 900 K. A 
number of bulk phosphide compounds can be produced, de- 
pending on the stoichiometric proportions of M/P used, but 
the Fe2P, CoP, and Ni2P products are the preferred phases 
under the experimental conditions used in this study. 
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FIG. 6. Temperature-programmed reduction of bulk nickel phos- 
phates at fi = 1 K/min (0.01667 Its' 1 ) (Ni-to-P ratios used in the prepa- 
ration are indicated). 



Properties and Preparation of Supported Phosphides 

In order to study the catalytic properties of the Fe, Co, 
and Ni samples high-surface-area materials were desired, 
and therefore, the phosphides were prepared in supported 
form using silica as the carrier. The preparation involved 
the same steps used in the synthesis of the bulk compounds. 
First, phosphates were prepared on the support and then 
were reduced by TPR. The final temperatures for prepar- 
ing hydrotreating catalysts were determined by noting the 
point where the intensity of the water signal returned to 
baseline in preliminary TPR test measurements on small 
samples. The final temperature for supported iron, cobalt, 
and nickel samples were 1000, 900, and 850 K, respectively 
(Fig. 8). Silica was selected as the carrier because in dehy- 
drated form it has few acid and base sites and is likely to 
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TABLE 3 
Characterization Results of Samples 





BET surface 


CO 






Metal site 




area 


uptake 




D c 


concentration 


Samples 


(S g /m 2 g" 1 ) 


Oxmol g* 1 ) 


(nm) 


(nm) 


(fimo\ g _1 ) 


Fe 2 P 


2.7 


3 


330 


38 


45 


j. i von * v7i loiiiwa 


97 


16 




23 


90 


Spent Fe 2 P/silica 


83 


0 




23 


90 


CoP 


3.1 


3 


310 


40 


51 


Fresh CoP/silica 


87 


16 




21 


72 


Spent CoP/silica 


90 


4 




21 


72 


Ni 2 P 


3.3 


4 


260 


36 


55 


Fresh Ni2P/siiica 


98 


15 




20 


67 


Spent Ni 2 P/silica 


100 


13 




20 


67 



offer minimal support interaction to affect the properties 
of the phosphides. Thus, the intrinsic activity of the phos- 
phides could be determined. 

Comparison of the TPR results obtained for the bulk 
samples and the corresponding supported samples (Fig. 8) 
provides strong evidence that the transformation of the 
phosphates to phosphides proceeded in the same manner 
for both bulk and supported samples. Aside from a slightly 
more pronounced low-temperature dehydration feature 
for the supported samples, the TPR traces for the sam- 
ples of iron and cobalt were very similar to those of the 
unsupported forms. The reduction trace of the supported 
nickel sample was somewhat simpler than that of the bulk 
sample, but occurred at essentially the same temperature 
range. 

The characterization results for the bulk and supported 
samples are reported in Table 3. The specific surface areas 
(S g ) of the bulk materials were low, approximately 3 m 2 g" 1 , 
while those of the supported materials were close to that 
of the support (90 m 2 g" 1 ). The experimental CO uptakes 
of the samples are reported in the third column of Table 3. 
They were low for the bulk materials but increased for the 
supported samples. 

The particle diameter of the bulk materials was calcu- 
lated by the equation Dp = 6/pS $ using their surface area 
and the bulk density (Table 2). The particle sizes for the 
bulk Fe2P, CoP, and Ni2P samples were similar, about 260- 
330 nm (Table 3). The crystallite sizes (D c ) were obtained 
from the Scherrer equation presented in the experimental 
section. The crystallite sizes for the bulk phosphides were 
again similar, about 36-40 nm. For all the bulk samples, 
Dp > D c , and this could be due to strain and disorder in 
the crystallites or crystallite agglomeration. The latter is 
probably the larger contributor, as no support was used to 
stabilize the samples. 

The last column in Table 3 reports the theoretical metal 
site concentration assuming that the samples were com- 
posed of uniform spherical particles. It was calculated from 



the equation 

Metal site concentration = 5 g • h • /, 

where S g is specific surface area, h is the surface metal 
atom density, and / is the fractional weight loading (e.g., 
grams of Fe 2 P/gram of catalyst) of the sample (Table 1). 
For the bulk samples the actual 5 g was employed, while for 
the supported samples it was calculated from the crystallite 
size using the equation S g = 6/pD c . The last factor, /, ac- 
counts for the loading of the active phase on the supported 
samples. 

In all cases the experimental CO uptake was consid- 
erably smaller (average 6.6% for the bulk samples; aver- 
age 20.5% for the fresh supported samples) than the the- 
oretically expected metal site concentration (Table 3) for 
a clean surface. This indicates that possibly the surface is 
blocked by some species that prevents adsorption, such 
as phosphorus or unreduced oxygen. In the case of tran- 
sition metal carbides it is found that oxygen uptakes are 
considerably higher than CO uptakes (42, 43), and this 
suggests that oxygen should be tried as a chemisorption 
probe. 

The effect of the heating rate on the peak temperature 
(T p ) associated with the reduction of supported catalysts 
(Fig. 9) was briefly examined in this work. The peak posi- 
tions for water formation shifted to a temperature about 
40-60 K higher as the heating rate (p) was increased from 
0.0167 (1 K min" 1 ) to 0.0833 K s" 1 (5 K min^ 1 ). According 
to temperature-programmed reaction theory (44), the peak 
temperature (T p ) is related to the heating rate (0) and the 
apparent activation energy (£ a ) (Table 4) by the Redhead 
equation (21nT p - ln/3 = E a //?r p + Constant) (44). The 
activation energies found for Fe2P/silica, CoP/silica, and 
Ni 2 P/silica were 200, 220, and 150 kJ mol" 1 , respectively. 
The results here are comparable with the activation energy 
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Apparent Activation Energy (£ a ) of Synthesis 
of the Supported Catalysts 



Samples 


£ a (kJ mol" 1 ) 


Fresh 14 wt% Fe 2 P/siIica 


200 


Fresh 9.4 wt% CoP/silica 


220 


Fresh 9.4 wt% Ni 2 P/siIica 


150 



of oxygen diffusion in the corresponding metal oxides, 
FeO, CoO, and NiO, which are 126, 144, and 166 kJmol" 1 , 
respectively, with preexponential factors of 1.4 x 10" 2 , 
2.15 x 10~ 3 , and 2 x 10~ 4 cm 2 s" 1 (45). The correspondence 
is reasonable, as many solid state transformations are gov- 
erned by diffusion processes (46). 



Catalytic Activity in Hydroprocessing 

Figures 10 and 11 present the HDS and HDN activities 
for the reactions of dibenzothiophene and quinoline, re- 
spectively. All three samples have high initial activities for 
HDS and HDN. However, except for the nickel sample in 
HDS, all catalysts undergo deactivation. The Fe2P/Si02 lost 
all HDS and HDN activity by 60 h, while the CoP/Si0 2 ap- 
peared to reach a baseline of about 32% HDS and 31% 
HDN at around 100 h. Only the Ni 2 P/Si0 2 had good, stable 
activity in HDS. The HDS sequence for the three samples 
was Ni 2 P/Si0 2 > CoP/Si0 2 > Fe 2 P/Si0 2 , while the HDN 
sequence was CoP/Si0 2 > Ni 2 P/Si0 2 > Fe 2 P/Si0 2 . Com- 
pared with a commercial Ni-Mo-S/y-Al 2 03 catalyst at the 
same conditions (24), Ni 2 P/Si0 2 had a higher HDS activity, 
with 90 versus 76% conversion, but a lower HDN activity, 
with 14 versus 38% conversion. The measurements were 
made on the basis of equal chemisorption sites loaded in 
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FIG. 10. Hydrodesulfurization performance of supported catalysts. 
(Basis: 35 /xmol of chemisorption sites). 
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FIG. 11. Hydrodenitrogenation performance of supported catalysts. 
(Basis: 35 /xmol of chemisorption sites). 



the reactor (35 /xmol for the phosphides, 33 /xmol for the 
sulfide). For the phosphides CO at room temperature was 
used for the chemisorption and for the sulfide O2 at dry 
ice/acetone temperature was used. It may be that the CO 
chemisorption underestimates sites on the phosphide, so 
the conversions should be taken just as an approximation 
of intrinsic activity. 

A listing of conversions and selectivities for all the re- 
actions is provided in Table 5. For dibenzothiophene the 
only product observed was biphenyL However, for quino- 
line a number of species were obtained, and these were 
categorized as HDN products and hydrogenation (HYD) 
products. For benzofuran, again only one product was ob- 
tained, ethylbenzene. The amount of benzofuran used was 
small (500 ppm) and the reaction is facile so is unlikely to 
interfere with the HDS and HDN reactions, as known from 
other studies with carbides, nitrides, and sulfides (47, 48). 
In the case of tetralin, at the reaction conditions the major 
species obtained was the dehydrogenation product naph- 
thalene, although small amounts of cis- and tVa/w-decalin 
were also observed. 

Turnover rates and areal rates measured at the baseline 
levels are listed in Table 6. The basis of comparison is now 
70 /xmol of sites, to make the numbers directly comparable 
to others published elsewhere (49). The procedure for 
conversion assumes a first-order reaction and is given in 
the footnotes of the table. The turnover rates are based on 
the experimentally determined CO uptakes on the fresh 
catalysts, while the areal rates are based on the calculated 
surface area of the phosphide crystallites using 5 g = 6/pD c 
and the weight loading factor /. The Fe 2 P catalyst deac- 
tivated completely and its rate is reported as zero. The 
CoP/Si0 2 had areal rates of 1.4 x 10 15 mol m~ 2 s~ l in HDS 
and2.1 x 10 15 mol nr 2 s' 1 in HDN, while the Ni 2 P/Si0 2 had 
rates of 3.4 x 10 15 mol nT 2 s" 1 in HDS and 1.4 x 10 15 mol 
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TABLES 

Product Distribution in Hydroprocessing" 



Conversion (%) Selectivity (%) 



Reactants 


Type 


Ni 2 P/Si0 2 


C0P/S1O2 


Product 


Ni 2 P/Si0 2 


CoP/Si0 2 


Dibenzothiophene 


HDS 


90 


32 


Biphenyl 


100 


100 


Quinoline 


HDN 


14 


31 


Propylcyclohexane 


4 


16 




HYD 


44 


49 


Propylbenzene 


9 


9 










5,6,7,8-Tetrahydroquinoline 


33 


30 










Orthopropylaniline 


31 


23 










1 ,2,3,4-Tetrahydroquinoline 


24 


22 


Benzofuran 


HDO 


35 


11 


Ethylbenzene 


100 


100 


Tetralin 


deHYD 


31 


6.5 


Naphthalene 


98.2 


91.5 




HYD 


0.6 


0.6 


frarts-Decalin 


0.6 


3.6 










cw-Decalin 


1.2 


4.9 



' Catalyst loaded was equivalent to 35 /xmol of CO uptake sites. 

TABLE 6 

Rates of HDS and HDN of the Supported Catalysts 



Turnover rate 0 



Areal rate* 



Specific rate 4 ' 



Volumetric rate* 



Sample 


(io- 3 


s" 1 ) 


(10 15 mol m" 


" 2 s- 1 ) 


(10 8 mol g" 


s- 1 ) 


(10 9 mol cm' 


-*s-i) 


HDS 


HDN 


HDS 


HDN 


HDS 


HDN 


HDS 


HDN 


Fe 2 P/Si0 2 


0 


0 


0 


0 


0 


0 


0 


0 


CoP/Si0 2 


0.81 


1.2 


1.4 


2.1 


1.3 


1.9 


4.8 


7.0 


Ni 2 P/Si0 2 


1.5 


0.60 


3.4 


1.4 


2.2 


0.9 


8.1 


0.33 



0 Calculated from r t = QX/S y where Q is the molar rate of reactant, X is the conversion, and S is the mole of sites loaded. The basis 
used was 70 ^mol. Conversions were adjusted using the first-order formula X 2 = 1 -(1 -X\ ) 5z/51 , where S 2 /S\ = 70 /xmol/35 fimol = 2, 
and gave higher X than reported in Table 5 (CoP: HDS, 54%, HDN, 52%; Ni 2 P: HDS, 99%, HDN, 26%). 

b Calculated from r A = QXN A /WS g f> where N A is Avogadro's number, W is the weight of catalyst, S g is the metal phosphide 
surface area, and / is the fractional loading of phosphide. 

c Calculated from r s =r t (CO uptake). 

d Calculated from r v =r s p, where p is the apparent density of the catalysts, ~0.37 g cm" 3 . 
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FIG. 12. X-ray diffraction patterns of iron phosphides and references. 
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FIG. 13. X-ray diffraction patterns of cobalt phosphides and refer- 



(a) Blank sample, Si0 2 ; (b) spent sample, Fe 2 P/Si0 2 ; (c) fresh sample, ences. (a) Blank sample, Si0 2 ; (b) spent sample, CoP/Si0 2 ; (c) fresh sam- 
Fe 2 P/Si0 2 ; (d) bulk sample, Fe 2 P; and (e) PDF 33-670 Fe 2 P (Ref. 41). pie, CoP/Si0 2 ; (d) bulk sample, CoP; and (e) PDF 29-497 CoP (Ref. 41). 
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m~ 2 s" 1 in HDN. In a recent study of the HDN on un- 
supported samples Stinner et al. report areal rates of 4.8 x 
10 16 mol m~ 2 $' 1 on C02P and 5.8 x 10 16 mol m" 2 s" 1 on Ni 2 P 
(26). The higher HDN rates on those samples can be attri- 
buted to several factors. First, Stinner et al. studied 
only HDN and employed a more reactive substrate (o- 
propylaniline rather than quinoline), which probably led to 
higher rates. They also carried out their tests without using 
sulfur compounds in the feed. Sulfur compounds can lead 
to competitive adsorption and lower rates. Finally, Stinner 
etal used low-surface-area (1-3 m 2 g~ l ) unsupported mate- 
rials in a low-conversion regime, which would tend to result 
in higher rates. Nevertheless, despite the differences in ex- 
perimental conditions the rates are of similar orders of 
magnitude. 

The XRD patterns of the silica support, fresh phosphide 
samples, spent samples, and PDF references are compared 
in Figs. 12-14. These results show that silica did not in- 
fluence the phases of the phosphides formed, which were 
the same as those obtained in the bulk materials Fe2P, 
CoP, and Ni2P The XRD patterns for the spent iron and 
cobalt samples were unchanged from the corresponding 
fresh samples, which shows that the Fe 2 P phase (Fig. 12) 
and CoP phase (Fig. 13) are stable during the hydrotreat- 
ing reaction. For the spent nickel sample, one more peak 
was observed (Fig. 14). Comparing this pattern with those 
of other nickel phosphide compounds, it was found that 
the peak (*) matched one due to Ni^Ps. This indicated 
that possibly a part of the Ni 2 P phase transformed to a 
Nii 2 P5-like phase during hydrotreating. XPS measurement 
for this sample (Table 7) showed a P/Ni ratio close to 1 
for the fresh sample, which decreased to 0.77 for the spent 
catalyst. This decrease may be related to the drop in HDN 
activity for this sample. The XPS analysis surprisingly also 
indicated a lack of S or N on the spent sample. The mea- 
surements were repeated on samples treated in He (just 
to remove volatile compounds), but the results were the 
same. This is likely due to lack of sensitivity by XPS at 
the levels expected for surface sulfur species in these sup- 
ported samples. Already the signals for Ni and P were 
very low. 

The results of posthydrotreating characterization are also 
presented in Table 3. For the iron sample, the CO uptake 
decreased to zero after the hydrotreating. The loss of the 
CO uptake is likely associated with the dramatic deactiva- 
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FIG. 14. X-ray diffraction patterns of nickel phosphides and ref- 
erences (a) Blank sample-Si0 2 ; (b) spent sample Ni 2 P/SiC>2; (c) fresh 
sample Ni 2 P/Si0 2 ; (d) bulk sample Ni 2 P; and (e) PDF3-953 Ni 2 P (Ref. 41). 

tion of the iron sample (Figs. 10c and 11c). As to the spent 
cobalt sample, the CO uptake also decreased substantially, 
but not to zero, as in the case of the iron sample. This was in 
line with the less complete deactivation of the catalyst. In 
the case of nickel, again a decrease in CO uptake was ob- 
served, but the decrease was much smaller than those of the 
iron and cobalt samples. Overall, the trend appears to indi- 
cate that the decrease in catalytic performance is related to 
the loss of active sites. 

The reduction in the number of active sites is appar- 
ently caused by sulfidation of the catalysts. This is cur- 
rently being investigated by extended X-ray absorption 
spectroscopy (EXAPS) of the fresh and spent samples 
(49). The results indicate that the catalysts are phase 
pure and sulfidation is restricted to the surface of the 
phosphide crystallites, as X-ray diffraction analysis shows 
that the phosphide phase is retained. The surface sulfida- 
tion process appears to be more pronounced in iron and 
cobalt phosphide, which show considerable deactivation 
in HDS and HDN. The nickel phosphide appears to be 
more tolerant to sulfur and thus retains high HDS activity. 
This is entirely reproducible. Its HDN activity is lowered, 
probably because the latter reaction is structure sensitive 
(49), and the site requirements for nitrogen removal are 
disrupted. 



TABLE 7 

XPS Results for Supported Nickel Samples 



Fresh 9.4 wt% 


Element 


Cls 


Ols 


Nls 


S2p 


P2p 


Ni 2p 


Si 2p 


Ni 2 P/silica 


Concentration (%) 


16.73 


55.42 


0 


0 


0.80 


0.83 


26.22 


Spent 9.4 wt% 


Element 


Cls 


Ols 


Nls 


S2p 


P2p 


Ni 2p 


Si 2p 


Ni 2 P/silica 


Concentration (%) 


19.24 


53.16 


0 


0 


0.36 


0.47 


26.77 
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CONCLUSIONS 

Pure Fe2P, CoP, and M2P phases were successfully syn- 
thesized by means of temperature-programmed reduction 
of the corresponding phosphates. The silica-supported 
forms of these samples were also successfully prepared, 
with retention of the active phase and increased CO up- 
take and BET surface area. It was found that Fe 2 P/Si02 
had good initial activity for HDS of dibenzothiophene and 
HDN of quinoline, but that this activity fell to zero in about 
60 h. The C0P/S1O2 catalyst also deactivated but appeared 
to reach a stable baseline of 32% HDS and 31% HDN 
conversion. Only the Ni 2 P/Si0 2 had a stable and high con- 
version in HDS of 90%, although its HDN was low at 
14%. The deactivation in all cases was associated with a 
decrease in the number of surface metal sites as titrated by 
the chemisorption of CO. For the most promising catalyst, 
Ni2P/Si0 2 , HDS was higher than that of a commercial cata- 
lyst, Ni-Mo-S/y-AkOs, based on equal sites loaded in the 
reactor, as measured by CO uptake for the phosphide and 
low-temperature O2 chemisorption for the sulfide. The de- 
velopment and improvement of this group of phosphides in 
the hydroprocessing field is a promising area of research. 
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Silica-supported molybdenum phosphide, MoP/SiCh, and tung- 
sten phosphide, WP/SiC>2, were prepared and characterized for their 
catalytic activity in hydrodenitrogenation (HDN) and hydrodesul- 
furization (HDS). The silica-supported phosphides were tested in 
the simultaneous hydroprocessing of quinoline and dibenzothio- 
phene at 3.1 MPa and 643 K and were compared to a MoS2/SiC>2 
sample. The supported phosphides had superior HDN and lower 
HDS activity compared to the sulfide under these conditions. The 
conversion levels of MoP/Si0 2 were similar to those of bulk MoP, 
while the WP/SiC>2 was somewhat less active than its bulk WP 
counterpart. However, when compared on a CO uptake basis, the 
bulk materials had a clear performance advantage. Solid state nu- 
clear magnetic resonance measurements indicated 31 P shifts from 
85% H3PO4 of 255 ppm for WP and 214 ppm for MoP. Compar- 
ison of freshly prepared and passivated samples indicated that 
the supported phosphides were susceptible to surface oxidation 
on exposure to the atmosphere as a result of their small particle 

Size. © 2002 Elsevier Science (USA) 

Key Words: hydrodenitrogenation; hydrodesulfurization; silica; 
tungsten phosphide; molybdenum phosphide; quinoline; dibenzo- 
thiophene; 31 PNMR. 



INTRODUCTION 

Transition metal phosphides like MoP (1-3), WP (4), 
Ni2P, and C02P (5, 6) constitute a new class of catalysts 
that have recently been reported to be active for hydropro- 
cessing reactions. This is currently an important area, as 
there is a need to develop better catalysts to meet severe 
environmental restrictions on the content of sulfur in trans- 
portation fuels. Fundamental understanding of the differ- 
ences between compounds like MoP and M0S2 is needed 
in order to design improved materials. The objectives of 
the present investigation are to study the preparation of 
MoP and WP deposited on a silica support and to compare 
their performance to the unsupported catalysts. A further 

1 To whom correspondence should be addressed. E-mail: oyama@ 
vt.edu. 



objective is to determine the applicability of 31 P magic- 
angle spinning-nuclear magnetic resonance (M AS-NMR) 
spectroscopy for the study of the bulk and supported phos- 
phide samples. Except for one study in our group on WP 
(7), this technique has not been applied to transition metal 
phosphides and the spectra presented here are the first to 
be reported in the literature. The NMR measurements in- 
dicate that the silica-supported MoP and WP samples are 
susceptible to oxidation on exposure to the atmosphere. 
However, complementary X-ray diffraction measurements 
show that at reaction conditions the principal supported 
phases remain phosphides. 

The use of silica for hydroprocessing catalysts has not 
been widely studied because weak interactions with the ac- 
tive phase generally lead to low dispersions (8). It is gener- 
ally thought that the surface of the silica support contains 
coordinatively saturated bridging oxide species and surface 
hydroxyls but is essentially neutral in character (9). This is 
in contrast to the alumina surface which contains hydroxyl 
groups with amphoteric character (10) and which interacts 
strongly with metal (11, 12), sulfide (13), phosphate (14), 
and phosphide (15) adsorbents. It has been reported that 
silica-supported C0M0 catalysts prepared by impregna- 
tion exhibit low hydrotreating activity compared to AI2O3- 
supported catalysts (16), and this has been attributed to 
the low dispersion of Co and Mo after calcination. How- 
ever, it has been found that low loading (2-4 wt%) samples 
of MoS2/SiC>2 can perform better than their AI2O3 coun- 
terparts (17, 18). In fact, silica-supported Co sulfide having 
7.4% loading has been recently proposed as an inexpensive 
alternative to conventional catalysts (19). We have under- 
taken this study to investigate the suitability of S1O2 as a 
support for MoP and WP. 

EXPERIMENTAL 

Catalyst Preparation and Characterization 

Catalysts were prepared by incipient wetness impreg- 
nation of a silica support with aqueous metal phosphate 
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precursors, followed by temperature-programmed reduc- 
tion (TPR) in flowing hydrogen. The silica support (Cabot, 
Cabosil L90, 99.8%) consisted of 0.2- to 0.3-/xm aggre- 
gates of spherical particles ~30 nm in diameter prepared 
by vapor-phase hydrolysis of SiCU in a H 2 -0 2 flame 
(20). Stoichiometric amounts of ammonium paramolyb- 
date, (NH 4 )6Mo 7 024 • 4H 2 0 (Aldrich, 99%), or ammo- 
nium metatungstate (NH 4 )6Wi 2 0 3 9 • xH 2 0 (Aldrich, 90%) 
were combined with ammonium phosphate (NH 4 ) 2 HP04 
(Aldrich, 99%) in distilled water, and the solution was used 
to impregnate the silica support by the incipient wetness 
technique. The concentration of Mo, W, and P in the solu- 
tion was 0.73 mol/dm 3 . The moist paste was calcined in air 
at 773 K for 6 h, ground and mixed, then pressed, broken, 
and sieved between 650- and 1180-/xm screens to form pel- 
letized supported phosphate materials. The molar loadings 
of molybdenum and tungsten monophosphide were held 
constant at 1.16 mmol each of Mo, W, and P per g of silica, 
leading to 13 wt% MoP/Si0 2 and 20 wt% WP/Si0 2 . Sulfi- 
dation of samples containing 5.0 wt% M0O3 and 13 wt% 
M0O3 (0.37 and 1.16 mmol of M0O3 per g of silica) led to 
5.6 and 17 wt% MoS 2 /Si0 2 , respectively. 

The precursor phosphates were loaded in quartz u-tube 
reactors and were reduced to phosphides using linear tem- 
perature ramps (0.0833 K s~ l , 5.0 K min" 1 ) in flowing hydro- 
gen. Hydrogen flow was set at 650 /xmols -1 g -1 of starting 
material (e.g. 300 cm 3 (NTP) min" 1 H 2 flow for a 0.300-g 
sample) for all syntheses. The final temperature was 850 K 
for the MoP/Si0 2 and 1000 K for the WP/Si0 2 , with this 
temperature held for 2 h to ensure completion of reaction. 
Species in the effluent were identified periodically (2-5 s) by 
a mass spectrometer (Dycor/Ametek MA100) to follow the 
progress of the reaction. Water was the only gaseous reac- 
tion product observed; no P 2 {m/z = 62) or PH3 (m jz = 34) 
were detected. Occasionally, though, traces of a yellow in- 
soluble deposit were formed downstream in the reactor 
following the higher temperature preparation of WP/Si0 2 . 
After reduction, the phosphide was cooled to room temper- 
ature under 65 /xmol s~ l of helium, and typically was pas- 
sivated progressively by flowing 78 /xmol s" 1 of 0.1 mol% 
0 2 /He overnight, followed by 13 /xmol s~ l of 0.5 mol% 
0 2 /He for 2 h. The reactor was then disconnected from the 
flow system, and ambient air was allowed to diffuse in for 
24 h before collecting the sample from the tube. 

Unsupported MoP and WP reference samples were also 
synthesized from phosphate precursors as reported earlier 
(1, 4). The bulk MoP was prepared by TPR to 913 K at 
0.0833 K s' 1 (5 K min" 1 ), and the bulk WP by TPR to 973 K 
at the same heating rate. 

MoS 2 /Si0 2 samples were prepared by sulfiding 5 wt% 
Mo0 3 /Si0 2 and 13 wt% M0O3/S1O2 for 2 h in a stream of 
10 mol% H 2 S/H 2 at 678 K, following a temperature ramp 
of 0.0833 K s" 1 (5 K min" 1 ). The sulfide sample used for hy- 
droprocessing was sulfided directly from the oxide within 



257 

the catalytic reactor at a flow rate of 98 /xmol s" 1 (150 cm 3 of 
NTP min" 1 ) of 10% H 2 S/H 2 at atmospheric pressure. The 
characterization of the sulfide samples by 0 2 chemisorp- 
tion also directly followed sulfidation of the oxide. For the 
chemisorption measurements on MoS 2 /Si0 2 , 0.3 g of oxide 
was sulfided in 195 /xmol s _1 of 10% H 2 S/H 2 employing the 
same temperature program used for the catalyst samples. 

All gases were supplied by Air Products. Hydrogen, he- 
lium, and carbon monoxide were all 99.999% pure and were 
passed through a water trap (Alltech) before coming into 
contact with the sample. Specialty mixtures, including 0.5% 
O z /He, 30% N 2 /He, and 10% H 2 S/H 2 , were used as re- 
ceived. 

The total dynamic CO uptakes of the phosphides were 
measured by passing pulses of CO (5.6 /xmol) in a stream 
of 40 /xmol s" 1 of helium over 0.3-g samples at room tem- 
perature. Passivated samples were characterized after re- 
reduction at 723 K in H 2 . The phosphides used in all our 
studies were used within a month of preparation, as it was 
found that longer term storage in closed vials resulted in 
a slow degradation, probably due to oxidation. Uptakes of 
0 2 for the sulfides were measured in an identical manner at 
dry ice/acetone temperature (195 K). The 0 2 uptakes were 
multiplied by a factor of 2 to account for dissociative ad- 
sorption on the surface. Surface area measurements were 
performed on air-exposed samples by N 2 adsorption using 
a Micromeritics 2000 sorption unit. 

X-ray diffraction (XRD) spectra were collected with a 
Scintag XDS-2000 X-ray diffractometer using Ni-filtered 
CuKa (k = 0.1541 nm) radiation. Scans were collected at a 
rate of 0.035° 20 s~ l and increments of 0.03° 20. Crystallite 
sizes were estimated from the diffraction patterns using the 
Scherrer equation, D h ki = KX/ficosO. The constant K was 
taken here to be 0.9, and the average of the strongest two 
reflections was used to determine the reported values (21). 

Nuclear Magnetic Resonance 

Characterization of the samples by 31 P MAS-NMR was 
carried out using a Bruker (MSL 300) instrument operat- 
ing at 300 MHz. A single-pulse method was used with a 
pulse delay of 3 /xs, an acquisition delay of 3 ms, a recycle 
delay of 1 s, and approximately 14,500 acquisitions. Mea- 
surements were taken at different sample rotation (SR) 
speeds. Chemical shifts are referenced to 85% H3PO4. Sam- 
ples were loaded into Kel-F capped zirconia rotors without 
exposure to the atmosphere in a glove box {in situ), or were 
passivated with 0.5% 0 2 /He and exposed to the atmosphere 
(ex situ). 

Reactivity Testing 

Hydroprocessing activity was measured by hydrotreat- 
ment of a model liquid containing quinoline and diben- 
zothiophene at 643 K (370°C) and 3.1 MPa (450 psig) in 
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Compositional Characteristics of Silica-Supported Catalysts 
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TABLE 1 
Composition of Model Oil 



N, S, O Delivery rate 
Component mol% wt% (wppm) (/xmolrr 1 ) 



Quinoline 2.6 1.9 2000 580 

Dibenzothiophene 1.7 1.7 3000 380 

Benzofuran 0.6 0.4 500 130 

0 6,100 





mmbl of metal g 1 


Metal atom nm" 2 


Sample* 


of Si0 2 


of Si0 2 


5.6% MoS 2 /Si0 2 


0.37 


2.44 


13% MoP/Si0 2 


1.16 


7.67 


20% WP/Si0 2 


1.16 


7.67 



Tetradecane 



68 



76 



15,500 



an upward-flow, fixed-bed reactor as described earlier (22). 
Compositional characteristics and delivery rates of the feed 
liquid are summarized in Table 1. With the exception of 
tetradecane of purity 99 mol% , which was purchased from 
Fisher, chemicals were obtained from Aldrich and had a 
purity level of 97+%. 

To start a reaction, samples were placed in the catalytic 
reactor and pretreated before starting the oil delivery. The 
phosphide catalysts were pretreated in 100 fimol s" 1 of hy- 
drogen at 723 K and 1-atm pressure for 2 h. The sulfide cata- 
lyst was pretreated by sulfiding the oxide directly in the cata- 
lytic reactor, with a 10 mol% H 2 S/H 2 mixture at 678 K and 
1-atm pressure for 2 h. After the pretreatment, the reactor 
conditions were set to 3.1 MPa (450 psig), a hydrogen flow 
rate of 100 Atmol s" 1 (150 cm 3 NTP min" 1 ) , a liquid feed rate 
of 0.0833 cm 3 s" 1 (5 cm 3 min" 1 ), and 643 K. The amounts of 
catalysts used in this study were 2.6 g of MoS2/Si0 2 , 1.4 g 
of MoP/Si0 2 , and 3.1 g of WP/Si0 2 . These amounts corre- 
sponded to roughly similar quantities (70, 70, and 58 /xmol, 
respectively) of chemisorption sites. For the Si0 2 blank, an 
arbitrary amount of 3 g was employed, as the material had 
negligible chemisorption of CO. The compositions of hy- 
droprocessed liquid samples were identified with a Hewlett 
Packard 5890A gas chromatograph (equipped with a CP 
Sil 5B column) on samples collected at 2- to 3-h intervals. 
Reaction products were identified through matching of re- 
tention times with commercially available chemicals. Char- 
acterization of catalysts after reaction by XRD was made 
after washing the collected samples in /z-hexane, drying in 
air, and mounting ~l-mm-thick powder films on top of a 
glass slide. 

RESULTS 

Preparation and Characterization 

Silica-supported molybdenum phosphide, tungsten phos- 
phide, and molybdenum sulfide catalysts were prepared by 
incipient wetness impregnation of a silica support with am- 
monium salts of Mo, W, and P dissolved in distilled water. 
Compositional characteristics of each sample are summa- 
rized in Table 2. 

Oxidic catalyst precursors were reduced to phosphide 
form by temperature-programmed reduction in hydro- 



a Weight percent MoP, WP, or MoS 2 . 

gen and were subsequently passivated. Prior to the re- 
activity measurements and ex situ CO uptake char- 
acterization, passivated phosphides were pretreated in 
hydrogen at 723 K for 2 h to remove the surface oxide 
layer. 

The temperature-programmed reduction (TPR) trace at 
5 K min- 1 of MoP/Si0 2 is compared to that of bulk MoP in 
Fig. 1. The predominant reduction features appeared at 680 
and 820 K for the supported sample and 770 and 920 K for 
the unsupported reference. X-ray diffraction (XRD) pat- 
terns of the fresh and spent MoP/Si0 2 samples are shown 
in Fig. 2 along with that of a bulk MoP reference. An in- 
sert depicts the crystal structure of MoP. The diffraction 
patterns of the supported catalysts showed only features 
of bulk MoP and amorphous silica (the broad peak near 
22° 20) both before and after reaction. 

The TPR traces of WP/Si0 2 and WP are shown in Fig. 3. 
Here, the predominant reduction peaks occurred at 890 
and 970 K for the supported material and 1010 K for the 
unsupported reference. XRD patterns of the fresh and 
spent WP/Si0 2 samples are shown in Fig. 4 along with 
that of a bulk WP reference. Again the insert depicts the 
crystal structure. The diffraction patterns of the supported 
catalyst showed features only of WP and amorphous silica 
both before and after reaction in hydrotreating. 
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FIG. 1. Temperature-programmed reduction profiles in the prepara- 
tion of MoP/Si0 2 and MoP at 0.0833 K s" 1 (5 K min' 1 ). 
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FIG. 2. X-ray diffraction spectra for fresh and spent MoP/Si02 
samples. 
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FIG. 4. X-ray diffraction spectra for fresh and spent WP/Si0 2 samples. 



TABLE 3 

Physical Characteristics of Silica-Supported Catalysts 





Surface Crystallite 


CO uptake* 


Surface 


Dis- 




area* 


size a 






metal 


persion* 


Sample 


(mV 1 ) 


(nm) 


Fresh 


Spent 


sites 1 * 


(%) 


Si0 2 blank 


91 




0 


0 






5.6% MoS 2 /Si0 2 


93 




27' 






5.2 


17% MoS 2 /Si0 2 






V 






0.7 


13% MoP/Si0 2 


50 


15/ 


50 


26 


110 


4.9 


20% WP/Si0 2 


62 


16* 


19 


6.5 


94 


2.0 


MoP 


8 


18/ 


15 


4.5 


130 


0.2 


WP 


10 


23* 


10 


11 


140 


0.2 



° Passivated catalysts. 

h CO uptake at 298 K in //mol g" 1 ; sample pretreated at 723 K in H 2 . 
c Calculated in /xmol g" 1 from formula in text. 
d CO uptake (fresh)/moI of metal. 

e 0 2 uptake at 78 K; sample presulfided at 678 K in 10% H 2 S/H 2 . 

* Average of (100) and (101) directions. 

* Average of (101 ) and (002) directions. 

The chemisorption and surface area properties of the 
catalysts are summarized in Table 3. The incorporation of 
phosphides onto the silica reduced the surface areas of the 
samples in both cases. The surface area of the sulfide was 
unaffected, presumably because of its low metal loading 
and low preparation temperature. The experimental CO 
uptakes dropped during the catalytic test. 

Crystallite sizes, theoretical CO uptakes, and dispersions 
were calculated for MoP and WP catalysts based on their 
crystallite size, density, and loading. These values are also 
reported in Table 3. The calculation of crystallite sizes (Dhki) 
with the Scherrer equation was described earlier, in the 
experimental section. Table 3 also reports the theoretical 
number of surface metal sites (/xmol g" 1 ), assuming that 
the samples are composed of uniform spherical particles. It 
was calculated from 



Surface sites = S g - h 



[i] 



where S % is the specific surface area (m 2 g _1 ) either mea- 
sured (for the bulk samples) or calculated (for the sup- 
ported samples) from the crystallite size and density by the 
equation 

S, = — — . [2] 

The quantity n is the areal density, the number of metal 
atoms per unit surface, and / is the fractional weight load- 
ing (g of MoP/g of catalyst or g of WP/g of catalyst) of the 
sample. The areal density is taken here to be 16.1 fxmol m~ 2 
for MoP, as calculated from the value of the surface area 
per metal atom of 10.3 A 2 , as reported previously (2). The 
corresponding values for WP were derived^ from its crys- 
tal structure to be 14.5 ^mol m~ 2 and 11.6 A 2 atom -1 . The 
density, p, was taken to be 7.50 g cm" 3 for MoP and 12.4 g 
cm~ 3 for WP (23). The crystallite sizes of the phosphides 
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supported on silica were smaller than the crystallite sizes of 
the bulk material. 

Dispersions were calculated from the fresh chemisorp- 
tion uptakes, assuming a stoichiometry of 1 CO per surface 
metal atom. These values are also reported in Table 3. The 
dispersions were generally low, indicating poor distribution 
of the components on the support. This is typical of the be- 
havior of metallic compounds on silica supports. 

Nuclear Magnetic Resonance Spectroscopy 

Characterization of the catalysts by 31 P MAS-NMR was 
carried out on both fresh (in situ) and passivated (ex situ) 
samples. Figure 5 shows the NMR spectra of the bulk 
WP at different sample rotation (SR) speeds. The in situ 
spectra (Fig. 5a) displayed a complex pattern in which the 
positions of only two peaks, at 255 and -6 ppm (arrows), 
were invariant with rotation. The ex situ spectra (Fig. 5b) 
again showed a number of peaks, with the same ones re- 
maining invariant, but with the intensity pattern changed 
considerably, and with the peak at low field becoming dom- 
inant and slightly shifted to -3 ppm. Figure 6 shows simi- 
lar NMR results for the WP/Si0 2 sample. Again, the in 
situ sample (Fig. 6a) had two peaks that were invariant 
with rotation at 254 and -5 ppm (arrows), while the ex 
situ sample (Fig. 6b) had an altered pattern with the peak 
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FIG. 5. 31 P MAS-NMR spectra: (a) WP in situ, (b) WP ex situ. 



at low field becoming dominant and slightly shifted to —3 
ppm. Table 4 gives a summary of the observed chemical 
shifts. 

Figure 7 shows the NMR spectra of the unsupported 
MoP sample at in situ and ex situ conditions. The in situ 
spectrum (Fig. 7a) displayed a complex pattern that was 
completely different from that of the bulk WP sample, but 
which again had two peaks, at 214 and -5 ppm (arrows), 
that did not shift with sample rotation. The ex situ sample 
(Fig. 7b) showed a similar pattern but the signals were much 
noisier and were overlaid on a sinusoidal baseline. Again, 
two invariant peaks were observed, with the low-field peak 
enhanced and shifted to -3 ppm. Figure 8 shows the re- 
sults for the MoP/Si0 2 sample. The in situ spectra (Fig. 8a) 
displayed considerable noise and an uneven baseline, with 
weak and unresolved features at 199 and -36 ppm (arrows). 
The ex situ spectra (Fig. 8b) were completely dominated by 
the low-field feature at -5 ppm, with the high-field peaks 
barely discernible. 

Reactivity 

The progress of the hydroprocessing reactions of diben- 
zothiophene and quinoline at 643 K and 3.1 MPa using 
MoP/Si0 2 and WP/Si0 2 are shown in Fig. 9. Both materials 
reached stable conversion levels after only 20 h of reaction. 
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DISCUSSION 

reparation and Characterization 

The reduction profile obtained in the preparation of the 
oP/Si0 2 catalyst is presented in Fig. 1 and is compared 

the reduction profile of bulk MoP taken under similar 
•nditions. The reduction trace of MoP/Si0 2 showed a fea- 
re at low temperature (370 K) and two peaks at high 
mperatures (680 and 820 K). The low-temperature fea- 
re was due to the dehydration of the sample, while the 
>h-temperature peaks were indicative of the reduction 
•ges of the phosphate. XRD analysis of the intermedi- 
i stage indicated that the first high-temperature peak, at 

K > was due t0 initial stages of reduction of Mo to a 
ver oxidation state, while the second peak, at 820 K was 
stated with the further reduction of both phosphorus 
i molybdenum to form MoP. These results will be pre- 
ited in a separate, more detailed study of the preparation 
). It can be seen that the synthesis reaction was substan- 
ly over at 850 K and this was used as the final reduction 
iperature for the preparation of the MoP/Si0 2 catalyst 
loldmg time of 2 h was used at this temperature to ensure 
apletion of the reduction. 

Tie bulk MoP reduction profile was obtained at the same 
tmg rate of 0.0833 K s" 1 (5 K min" 1 ) as the supported 
tple and also showed two high-temperature reduction 
cesses. However, the first peak was unresolved and both 
ks were shifted to higher temperatures (770 and 920 K) 
» was likely due to the presence of larger particles react- 
in a shell-to-core manner and inhibition of reaction by 
it formed by the reaction (25). 

he TPR traces (Fig. 1) of the molybdenum phosphide 
pies can be analyzed to obtain the average rate of for- 
lon of the product phosphide, and hence the rate of 
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consumption of hydrogen. From the temperature program 
and duration of the main reduction peaks (full width at 
half maximum = 65 K) it is possible to extract the time of 
reaction (800 s). This information, combined with knowl- 
edge of the total amount of sample in the preparation, 
yields the average rate of hydrogen consumption. Assum- 
ing that the stoichiometry of the supported material is 
similar to that of the bulk, as determined gravimetrically, 
M0PO5.5 + 5.5H 2 -> MoP + 5.5H 2 0, the rate of hydrogen 
consumption for the supported sample can be calculated 
to be 6.9 izmol g" 1 s" 1 . This compares to the flow rate of 
hydrogen of 650 /zmol g~ 3 s _l , which is a factor of 94 times 
larger. Therefore, the reaction is unlikely to be inhibited by 
the water produced by the reaction. 

For the unsupported sample (full width at half maxi- 
mum = 120 K, or 1500 s) the corresponding rate of hydrogen 
consumption is 28 ^mol g _1 s _1 . The flow rate of hydrogen 
is now only 24 times larger, and the possibility of product 
inhibition by the water formed in the reaction is larger. 

Another effect which is also likely to be contributing to 
the retardation of the reduction peaks in the bulk mate- 
rial is shell-to-core reduction. The average particle size of 
the bulk starting material can be calculated using the ex- 
pression diameter = 6/surface area • density (where surface 
area = 1.3 m 2 g"\p = 3.54 cm 3 g _1 ) to be 1.3 /xm. This was 
likely to be much larger than the particle size of the sup- 
ported phosphate precursor. Since the rate of reaction in 
the shrinking-core model depends on the square of the par- 
ticle size (26, 27), this probably also contributed to the ob- 
served shift in the temperature of reaction between the bulk 
and silica-supported materials. 

Assuming the stoichiometry of bulk reduction indicated 
above the final stoichiometry of the phosphides is MoP and 
WP. As is discussed below this is supported by the XRD 
and NMR results. These compounds have a metallic nature 
(28), so it is expected that the oxidation states of Mo and W 
will be close to zero. However, it is likely that some charge 
transfer will occur to the electronegative P atoms, so in fact 
the Mo and W should be slightly oxidized. 

The XRD spectra of the bulk MoP and supported MoP 
samples prepared at 850 K, before and after the hydrotreat- 
ing reactions, are shown in Fig. 2. The bulk MoP had a 
pattern which corresponded exactly to that reported in the 
literature (powder diffraction file 24-771). MoP has a hexa- 
gonal WC-type structure (strukturbericht designation, fl h ) 
with space group P~ bm2 and lattice parameters a 0 = 322 pm 
and c 0 = 319 pm (29). The insert shows a schematic of the 
structure, which is based on a stacked arrangement of tri- 
gonal prisms. The supported MoP/Si0 2 showed a broad fea- 
ture at —22° 20 produced by the amorphous silica support, 
and several distinct peaks which matched those of the bulk 
MoP phase. 

The reduction profile for the preparation of WP/Si0 2 
catalyst is presented in Fig. 3 and is compared to the reduc- 



tion profile for bulk WP collected under similar conditions. 
The reduction trace of WP/Si0 2 showed a small feature at 
370 K, a slight rise at 630 K, and a large peak at 890 K 
with a high-temperature shoulder at 970 K. The feature 
at 370 K was due to dehydration of the sample, while the 
high-temperature peaks were due to formation of phos- 
phide. The slight rise at 630 K may have been due to small 
quantities of highly dispersed species, while the shoulder 
at high temperature may have arisen from larger particles 
which reduced more like the bulk WP. It is clear from the 
figure that the reduction process was essentially completed 
at 1000 K, which was used as the maximum temperature in 
the preparation of the WP catalyst. Again, a hold time of 
2 h was used to ensure complete reduction of the sample. 

The reduction profile of the bulk WP also showed the 
slight dehydration feature at 370 K and a high-temperature 
reduction peak at 1010 K. Similarly to the case of supported 
and unsupported MoP, the largest peak was shifted to higher 
temperature, likely as a result of the inhibitive action of 
water vapor and the presence of larger particles in the bulk 
material. 

The XRD spectra of the bulk WP and supported WP 
samples prepared at 1000 K, before and after hydrotreat- 
ing reactions, are shown in Fig. 4. The bulk WP had a 
pattern corresponding exactly to that reported in the lit- 
erature (powder diffraction file 29-1364). WP adopts the 
orthorhombic MnP type structure (strukturbericht desig- 
nation, fl 3 i) with space group P nma and lattice parameters 
a 0 = 571 pm, b 0 = 325 pm, and c 0 = 623 pm (30). The 
insert shows a schematic of the structure, which like MoP 
is based on a stacked arrangement of trigonal prisms. How- 
ever, in WP the second row is displaced laterally half a lat- 
tice spacing and the phosphorus atoms are linked in chains. 
The supported WP/Si0 2 showed a broad feature at ~22° 20 
due to the silica support, and several distinct peaks which 
matched those of the WP phase. 

The ex situ chemisorption characteristics and BET sur- 
face areas of catalyst samples are summarized in Table 3. 
The surface area of the MoS 2 /Si0 2 catalyst (93 m 2 g" 1 ) pre- 
cursor was close to that of the support (91 m 2 g" 1 ), as its 
loading was low. The surface areas of the MoP/Si0 2 (50 m 2 
g" 1 ) and WP/Si0 2 (62 m 2 g" 1 ) samples we re significantly de- 
creased from that of the support. This is a common finding 
for phosphorus-promoted, supported catalysts (31, 32), and 
probably occurs as a result of solid state reactions or pore 
clogging promoted by phosphorus compounds. In the case 
of the phosphides studied here, the use of elevated temper- 
atures during sample reduction, 850 K for MoP/Si0 2 and 
1000 K for WP/Si0 2 , may have contributed to the decrease 
in surface area by sintering of the silica support. Chemisorp- 
tion quantities of 0 2 and CO for the fresh catalysts were 
moderate and were found to decrease somewhat after re- 
action for the phosphides. Measurements taken after re- 
action were carried out on recovered samples washed in 
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hexane and re-reduced at 723 K for 2 h. The CO uptakes of 
MoP/Si0 2 and WP/Si0 2 and 0 2 uptake of MoS 2 /Si0 2 were 
used to calculate dispersions based on stoichiometrics for 
adsorption of one CO molecule and one O atom per sur- 
face metal atom. The dispersions are reported in Table 3. 
They were low, 5.2% for 5.6% MoS 2 /Si0 2 , 0.7% for 17% 
MoS 2 /Si0 2 , 4.9% for MoP/Si0 2 , and 2.0% for WP/Si0 2 , 
consistent with the observation of crystalline phases for the 
phosphides. 

The very low dispersion of the 17% MoS 2 /Si0 2 was prob- 
ably due to the formation of large domains exposing in- 
active basal planes. For this reason the low-loading 5.6% 
MoS 2 /Si0 2 sample was used for the reactivity measure- 
ments. 

The crystallite sizes of the phosphide samples were deter- 
mined by Scherrer analysis and are reported in Table 3. As 
indicated in the table, the crystallite sizes were taken as av- 
erages over two principal directions. Although not shown, 
the crystallite dimensions in those directions were not sub- 
stantially different (<5%), indicating that the particle mor- 
phology was globular. 

The crystallite sizes can be used to calculate the theo- 
retical number of surface metal sites (fimo\ g _1 ) based on 
clean surfaces and spherical particles, and these are also re- 
ported in Table 3. For the unsupported materials the num- 
ber of surface sites were considerably larger than the CO 
uptakes, indicating that a large proportion of the surface 
(89% for MoP, 93% for WP) was unable to chemisorb CO. 
This may have been due to site blockage by phosphorus or 
strongly held oxygen. For the supported samples the num- 
ber of sites was closer to the actual uptakes, indicating that 
a smaller proportion of their surfaces was blocked (55% for 
MoP/Si0 2 , 80% for WP/Si0 2 ). 

Nuclear Magnetic Resonance 

The 3l P MAS-NMR spectra (Figs. 5-8) are the first re- 
ported for supported metallic phosphides in the literature. 
The spectrum of the bulk WP (Fig. 5a) at zero rotation 
showed an asymmetric envelope characteristic of consider- 
able electric-field anisotropy surrounding the phosphorus 
nucleus. This can be understood from the crystal struc- 
ture (Fig. 4), which shows that the phosphorus atoms re- 
side in sites which look different, when approached from 
the jc, y, or z directions. When the samples were rotated, 
the anisotropy was folded over to the spinning side bands, 
resulting in a complex intensity pattern. Nevertheless, crys- 
tallographically all the phosphorus atoms were equivalent, 
and the position of only one major peak, at 255 ppm, was 
invariant with rotation. A smaller peak at -6 ppm was due 
to a small amount of phosphate impurity. Chemical shifts 
in this range are typical for phosphates (14). For the passi- 
vated WP sample (Fig. 5b) this peak shifted to -3 ppm and 
became the dominant feature of the spectrum. The pas- 
sivation procedure for the NMR samples was less gentle 



TABLE 4 



Summary of Chemical Shifts for Phosphide Catalysts 0 



Sample 


Phosphide (ppm) 


Phosphate (ppm) 


In situ 


Passivated 


In situ 


Passivated 


WP 


255 


254 


-6 


-3 


WP/Si0 2 


254 


254 


-5 


-3 


MoP 


214 


214 


-5 


-3 


MoP/SiQ 2 


199 


None 


-36 


-5 



a Reference: 85% H3PO4. 



(0.5% 0 2 /He followed by air exposure) than for the cata- 
lysts (gradual increase in 0 2 concentration) and probably 
resulted in more oxidation. There is also the likelihood that 
the intensity of the insulator phosphate signal is intrinsi- 
cally stronger than that of the metallic phosphide, where 
the phosphorus is surrounded by conduction electrons. 

The spectra of the supported WP/Si0 2 samples (Fig. 6) 
show features very similar to the bulk references (Fig. 5). 
The in situ pattern (Fig. 6a) indicated a larger amount 
of phosphate than in the bulk reference. This was likely 
due to the smaller crystallite size in the supported ma- 
terials (Table 3), which resulted in easier oxidation. In 
passivated samples the phosphate peak dominated the pat- 
tern (Fig. 6b), but the phosphide signature was clearly 
present. XRD analysis of this sample (Fig. 4) indicated 
the presence of crystalline WP, supporting the contention 
that the phosphate NMR signal is intrinsically strong com- 
pared to the phosphide, and not a good indicator of the 
real amounts of each phase. All the chemical shifts are 
summarized in Table 4. There are small shifts in the lines, 
but all can be readily assigned to phosphate or phosphide 
compounds. 

The spectrum of the bulk MoP at zero rotation (Fig. 7a) 
also showed an asymmetric envelope, but different from 
that of WP. This is because MoP has a different structure 
(Fig. 2) than WP, but still has considerable anisotropy sur- 
rounding the phosphorus atoms. Rotation of the sample 
again produced a complex pattern of spinning side bands 
caused by the anisotropy, but one main peak, at 214 ppm, 
remained invariant. A second smaller peak, at -5 ppm, 
was due to the presence of some phosphate contaminant. 
Again, as with WP it is expected that the intensity of this 
phosphate peak is much stronger than that of the metal- 
lic phosphide signal, so the relative size of the peaks is 
not indicative of the amounts of each substance. When the 
sample was exposed to air (Fig. 7b) the phosphate peak 
shifted slightly, to -3 ppm, and grew in intensity. The noise 
level also went up substantially. For the case of the sup- 
ported MoP/Si0 2 at in situ conditions (Fig. 8a) the noise 
level was also large, and the peaks for the phosphide and 
phosphate were barely discernible. In the case of the pas- 
sivated sample (Fig. 8b) the phosphate signal completely 



264 



CLARK, WANG, AND OYAMA 



dominated the spectrum. It is likely that the MoP/Si0 2 is 
more easily oxidizable than the WP/Si02. Again, a sum- 
mary of chemical shifts is given in Table 4. In conclu- 
sion, for both WP and MoP the NMR spectra showed 
that supporting the materials on S1O2 allowed retention 
of the phosphide phases. Passivation of the samples re- 
sulted in some oxidation, with a greater severity for the MoP 
system. 

Reactivity 

The activity of the catalysts in the hydrodesulfurization 
(HDS) of dibenzothiophene and the hydrodenitrogenation 
(HDN) of quinoline was evaluated in a three-phase flow re- 
actor at 643 K and 3.1 MPa. Figure 9 shows the conversions 
as a function of time on the silica-supported MoP and WP 
in HDS and HDN. Both materials demonstrated good sta- 
bility in their activity levels. It should be noted that those 
results were obtained with 58 /Ltmol of sites of WP/S1O2 and 
70 /zmol of sites of MoP/Si0 2 loaded in the reactor. The 
conversions were corrected to equal sites assuming a first- 
order kinetic expression. Denoting X\ as the conversion 
corresponding to a quantity of sites Si, and X2 as the con- 
version corresponding to Si sites, then the quantities are 
related by 

ln(l-X 2 )= (J^j -ln(l-Xi). 

The results, corrected to 70 /xmol of sites, for the supported 
MoP/Si0 2 , WP/Si0 2 , and MoS 2 /Si0 2 are summarized in 
Table 5, along with results from bulk MoP and WP and 
the silica support. The conversion levels of MoP/Si0 2 were 
similar to those of the bulk MoP (1). However, adjustment 
by the CO uptake made the bulk appear better. This was 
a consequence of the low CO uptakes of the bulk material 
compared to the supported sample. Because of the low sur- 

TABLE 5 



Summary of Hydroprocessing of Model Feedstock 
at 643 K and 3.1 MPa 



Sample" 


HDN" (%) 


HYD' (%) 


HDS'' (%) 


Si0 2 


2 


32 


1 


5.6% MoS 2 /Si0 2 


11 


45 


31 


13% MoP/Si0 2 


31 


45 


20 


20% WP/SiCV 


46 


46 


16 


MoP r 


54 


34 


24 


WP 


88 


33 


92 



" Reactor catalyst loading was 70 /xmol of sites except for WP/Si0 2 
(58 /xmol) and bulk samples (30 m 2 ). 
h Hydrodenitrogenation of quinoline. 
'* Hydrogenation of quinoline to N-containing products. 
d Hydrodesulfurization of dibenzothiophene. 

r Percentage of HDN and HDS adjusted to 70 //mol of sites according 
to first-order equation given in text. 



face area of the bulk MoP, a strong possibility exists that a 
large fraction of the active surface is underestimated. The 
WP/Si02 sample differed most substantially from its un- 
supported counterpart (4) by having lower HDS activity. 
Again, this was due to the low CO uptake of the bulk ma- 
terial. Compared to the MoS 2 /Si0 2 both phosphides had 
significantly higher HDN, but lower HDS. 

The low apparent activity of the supported phosphides 
compared to the bulk references was surprising. It may 
have been that passivation resulted in overoxidation of the 
catalysts, as shown by NMR. Re-reduction of the catalysts 
produced the phosphides again, but perhaps in a different, 
less-active form from the bulk phosphides. This form may 
have had phosphides that were phosphorus deficient be- 
cause of segregation. Nevertheless, the samples were com- 
pletely stable to reaction conditions, as clearly shown by the 
XRD patterns after reaction (Figs. 2, 4). 

In recent work (3) the group of Prins studied the HDN 
on a series of bulk phosphides at conditions similar to those 
employed in this work (3.0 MPa and 643 K). In the presence 
of sulfur they found HDN conversions of 10% for MoP and 
30% for WP, which are somewhat lower than those found 
in this study, 31 and 46%, respectively. The comparison is 
not exact because the previous work used 0-propylaniline 
instead of quinoline, bulk catalysts instead of supported 
catalysts, and slightly different catalyst weights. Consider- 
ing these differences the results are comparable. 

An interesting conclusion of this study is that sample mor- 
phology is important, especially in the case of a weakly in- 
teracting support like silica. In the case of molybdenum 
sulfide the effect of loading showed that low concentra- 
tions were required to produce samples of sufficient active 
surface. This can be understood from the layered structure 
of M0S2, which has inactive basal planes and active edge 
planes. As the size of the crystallites grow with loading, 
they do so at the expense of the desirable side planes. How- 
ever, for molybdenum and tungsten phosphide, which do 
not have a layered structure, the smaller interaction with 
silica probably produces particles of a globular morpho- 
logy whose site density is much higher. This represents a 
special opportunity for these new materials, which could be 
exploited by the use of better supports than the silica used 
in this study. 

CONCLUSIONS 

MoP and WP were prepared on a silica support. They 
were moderately active catalysts for the hydrodenitrogena- 
tion of quinoline and the hydrodesulfurization of diben- 
zothiophene in a complex mixture. The silica-supported 
catalysts had similar temperature-programmed reaction 
profiles as those of the bulk samples, while having reduced 
hydroprocessing activities compared to the bulk materials. 
The smaller activity was partly due to the inaccuracy of 
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counting sites by CO chemisorption on the low-surface- 
area bulk materials. It could also have been the result of 
oxidation of the small particles of the phosphide on the 
support on exposure to air. A low loading was required to 
generate a MoS2/Si02 samples of sufficient O2 uptake for 
catalytic testing. 

The materials were characterized using 31 P magic-angle 
spinning-nuclear magnetic resonance spectroscopy. WP 
and MoP gave isotropic signals shifted 255 and 214 ppm 
from 85% H3PO4, respectively. The presence of oxide was 
indicated by signals between —3 and -6 ppm whose inten- 
sity increased substantially with passivation. 
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Abstract 

Alumina-supported molybdenum phosphide hydroprocessing catalysts, MoP/y-Al203, with weight loadings from 3.5 to 39 wt% were 
prepared by temperature-programmed reduction of alumina-supported molybdenum phosphate precursors. The precursors were obtained 
by incipient wetness impregnation of the support with aqueous molar solutions of ammonium paramolybdate and ammonium phosphate 
(Mo/P = 1/ 1 ). Effects of loading, reduction temperature, and heating rate on the catalysts were studied, and the samples were characterized 
by CO chemisorption, BET surface area, and X-ray diffraction (XRD) measurements. Compared to the bulk MoP and MoP/SiC>2 systems 
(which are similar to one another), the MoP/y-Al203 material showed different behavior. Whereas calcined bulk and silica-supported 
molybdenum phosphates form Mo-P amorphous glasses which reduce directly to MoP, for the alumina-supported phosphates at high loadings 
(e.g., 13 wt% MoP/y-A^C^), a M0O3 component is visible which reduces sequentially with temperature to M0O2 and Mo metal, and then 
transforms to MoP. The difference in behavior between these systems is attributed to the formation of a strongly bound amorphous phosphate 
surface layer on the alumina support. Alumina has a strong affinity for phosphates, and appears to initially abstract a substantial amount of 
phosphorus from the stoichiometric molybdenum phosphate mixture. The alumina then releases the phosphorus at high temperature, allowing 
it to recombine with Mo metal to form MoP and generating active sites. The MoP/y-A^C^ catalysts performed well in hydroprocessing of a 
simulated distillate containing dibenzothiophene and quinoline at conditions representative of industrial hydroprocessing (643 K, 3.1 MPa). 
For the 13 wt% sample the hydrodesulfurization conversion was 57% and hydrodenitrogenation conversion was 62%. Catalytic activity, 
based on equal chemisorption sites loaded in the reactor (70 umol), was generally independent of the amount of MoP deposited on the 
alumina surface, independent of the presence of X-ray visible molybdenum phases, and was associated with a relatively high temperature 
reduction peak found at all loading levels of M0P/AI2O3 but not in bulk MoP. We conclude that the high-temperature peak gives rise to 
highly dispersed MoP which is responsible for the bulk of the CO titration sites and the catalytic activity, and that the large MoP particles 
visible by XRD make smaller contributions. 
© 2003 Elsevier Inc. All rights reserved. 

Keywords: Molybdenum; Phosphide; Phosphorus; Alumina; Hydrodesulfurization; Hydrodenitrogenation; Dibenzothiophene; Quinoline 



1. Introduction 

Catalytic hydrotreating is an important refining process 
used to remove sulfur, nitrogen, and oxygen from hydrocar- 
bon fuel liquids. Research in our laboratory has identified 
molybdenum phosphide, MoP, as an active and stable hy- 
droprocessing catalyst [1—3]. Additional work by the group 
of Prins [4] has shown that MoP has six times higher ac- 
tivity in the HDN of orthopropylaniline than M0S2/AI2O3. 
Studies in the group of Bussell have shown that M0P/S1O2 
has four times higher activity in the HDS of thiophene than 
M0S2/AI2O3 [5]. In this paper we report the synthesis and 
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catalytic properties of MoP deposited on a y -alumina sup- 
port. 

MoP has a hexagonal WC-type structure, space group = 
P6m2, with lattice parameters a Q = 322 pm and c Q = 
3 1 9 pm [6]. MoP is a metallic conductor with properties sim- 
ilar to ordinary intermetallic compounds [7]. Fig. 1 shows 
the WC structure type adopted by MoP. The WC structure 
type, along with the NiAs and WP structures, characterize 
the majority of transition metal monophosphides. The WP 
structure type is similar to NiAs, but distortions in the (001) 
plane form a rectangular net, resulting in an orthorhom- 
bic structure [8]. A recent review summarizes the structure, 
preparation, and catalytic activity of phosphides [9]. 

Removal of heteroatoms from fuel liquids is important 
because sulfur and nitrogen compounds corrode engine com- 
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Table 1 

Compositional characteristics of M0P/AI2O3 samples 



• Mo • P 

Fig. 1. Tungsten carbide (WC) structure type adopted by molybdenum 
phosphide, MoP. 



ponents, pollute the air as SO* and NO*, cause gelling of 
untreated liquids, and foul downstream reforming catalysts. 
The problem is becoming increasingly important as gov- 
ernment regulations become stricter [10] and the quality of 
petroleum feedstocks declines [11]. Although hydroprocess- 
ing catalysts have been widely studied, their basic composi- 
tions have changed little for many decades. They are typ- 
ically alumina-supported molybdenum sulfides with nickel 
or cobalt promoters, NiMoS/Al 2 0 3 or C0M0S/AI2O3 [12]. 
The catalysts described in this work are completely differ- 
ent, consisting of a compound of Mo and P with metallic 
rather than insulating or semiconducting properties. 

In recent years, phosphorus has come to be used as a 
secondary promoter in commercial catalysts [13]. The use 
of phosphorus in hydroprocessing catalyst formulations was 
introduced in a series of patents dating back to the 1950s 
[14-16]. The phosphorus in these catalysts is dispersed on 
the catalyst surface in oxidized form, as shown by studies 
with X-ray photoelectron spectroscopy [17], nuclear mag- 
netic resonance [18], and infrared spectroscopy [19]. The 
effect of phosphorus depends on its concentration and phys- 
ical form, but in general it has a larger positive effect on the 
HDN performance of Ni-Mo sulfides than in the HDS per- 
formance of Co-Mo sulfide catalysts [20]. 



2. Experimental 

Catalysts were prepared by impregnation and reduction. 
Desired quantities of ammonium paramolybdate tetrahy- 
drate, (NH 4 )6Mo 7 0 2 4 ■ 4H 2 0 (Aldrich, 99%), and ammo- 
nium phosphate, (NH4) 2 HP0 4 (Aldrich, 99%), were dis- 
solved in 1 : 1 stoichiometric proportions of metal to phos- 
phorus to form clear solutions in distilled water. These so- 
lutions were used to fill the pores of a y -alumina support 
(Degussa, Aluminumoxid C) by incipient wetness impreg- 
nation. The moist paste was calcined in air at 773 K for 6 h, 
ground with a mortar and pestle, pressed, broken, and sieved 
to particles between 650 and 1 1 80 urn in diameter for use 



wt% MoP 


mmolMoPg 1 A1 2 C>3 


Molecules MoP nm 2 AI2O3 


3.5 


0.29 


1.9 


6.8 


0.58 


3.8 


13 


1.16 


7.7 


23 


2.3 


15 


39 


3.8 


25 



(mesh size 14/20). Samples were prepared with loadings be- 
tween 0 and 39 wt% MoP (0 and 3.8 mmol g" 1 A1 2 0 3 ). Spe- 
cific compositional characteristics for the series of catalysts 
are reported in Table 1 . The highest loading sample was pre- 
pared by double impregnation, due to solubility limitations 
of the metal and phosphate salts in distilled water. Much of 
this study used the 13 wt% (1.16 mmol MoP g" 1 Al 2 03) 
MoP/Al 2 03 material. Reference samples containing only 
1.16 mmol Mo g~ l A1 2 0 3 (10.0 wt% Mo) or 1.16 mmol 
P g~ l Al 2 03 (3.5 wt% P) were prepared in a similar man- 
ner to the catalysts, while aluminum phosphate, AIPO4, was 
purchased from a commercial source (Alfa-Aesar Division 
of Johnson-Matthey, 99.999%). 

Temperature-programmed reductions were carried out on 
pelletized catalyst samples (typically 0.3 g) placed in quartz 
u-tube reactors. The samples were heated with linear temper- 
ature ramps in flowing hydrogen (Airco, 99.999%) to reduce 
the metal phosphate to metal phosphide. Unless otherwise 
noted, the space velocity was set at 650 umol H 2 s~ l g -1 
of sample loaded (1000 cm 3 (NTP) min" 1 g" 1 ), and the 
heating rate was 0.0833 Ks*" 1 (5 Krnin" 1 ). The effluent 
was monitored continuously by a mass spectrometer (Ame- 
tek/Dycor, MA 100). All samples were flushed in hydrogen 
flow for 2 h before starting the temperature program, to 
dry the sample and to provide a starting baseline in the 
mass spectrometer. After reduction, the product was cooled 
to room temperature (RT) under helium (Airco, 99.999%) 
flow, and passivated progressively with 0.1% 0 2 /He, 0.5% 
0 2 /He (Airco), and slow diffusion of air to the sample 
for 24 h. 

The total dynamic CO uptake (Air Products, 99.999%) 
was measured with a mass spectrometer by passing 5.6 umol 
pulses over the sample in 65 umol s" 1 helium flow at RT un- 
til constant area peaks were obtained. Measurements made 
on freshly reduced samples are referred to here as in situ 
CO uptakes while measurements taken on passivated sam- 
ples rereduced in H 2 for 2 h at 723 K are referred to here as 
ex situ CO uptakes. BET specific surface areas (5 g ) of pas- 
sivated samples were measured with a Micromeritics 2000 
sorption unit. 

X-ray diffraction (XRD) patterns were obtained at a scan 
rate of 0.035° 2(9 s" 1 using a Scintag XDS-2000 diffrac- 
tometer with a Ni-filtered Cu-K a (X = 0. 1541 nm) radiation 
source. Crystallite sizes were estimated using the equation of 
Scherrer, D = Kk/ficosO. The constant K was taken here 
to be 0.9, k is the wavelength of radiation, p is the peak 
width in radians, and 0 is the angle of the reflection. The 
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average of the first three reflections in MoP was used to cal- 
culate an average particle diameter [21]. Activation energies 
for reduction reactions were estimated by the heating rate 
variation method of Redhead [22]. 

Hydrotreating activity was evaluated with a model distil- 
late containing 2000 wppm N as quinoline (Aldrich, 99%), 
500 wppm O as benzofuran (Aldrich, 99%), 3000 wppm S as 
dibenzothiophene (Aldrich, 99%), 20 wt% tetralin (Aldrich, 
99%), and balance tetradecane (Fisher, 99%). The catalytic 
test was carried out in an upflow, three-phase, fixed-bed re- 
actor as described earlier [23]. An amount of catalyst cor- 
responding to 70 umol of CO sites (measured ex situ) was 
loaded for each reaction. To start the reaction, catalysts were 
pretreated in flowing hydrogen (100 umol s" 1 ) at 723 K and 
1 atm pressure for 2 h. The pressure was then increased to 
3.1 MPa (450 psig), the temperature was dropped to 643 K 
(370 °C), and a liquid feed rate of 5 cm 3 min" 1 (delivering 
9.7 umol min -1 of quinoline and 6.3 umolmin" 1 of diben- 
zothiophene) was initiated. Analysis of hydroprocessing liq- 
uids was performed with a gas chromatograph (Hewlett- 
Packard, 5890A) on samples collected at 2- to 3-h intervals 
for a time period of 60-80 h, until a stable steady state was 
reached. In this report, hydrodesulfurization (HDS) is de- 
fined as the percentage disappearance of dibenzothiophene 
(as no sulfur-containing intermediates were observed), hy- 
drogenation (HYD) is defined as the percentage of original 
quinoline converted to hydrogenated quinolines, and hydro- 
denitrogenation (HDN) is defined as the percentage of nitro- 
gen removed from the quinoline delivered to the catalyst. 



3. Results and discussion 

3. 7. Synthesis of MoP /y -Al 2 0 3 

Calcined catalyst precursors were reduced to M0P/AI2O3 
catalysts by temperature-programmed reduction in hydro- 
gen at 0.0833 Ks~' (5 Kmin" 1 ). Corresponding reduction 
profiles for 0-23 wt% M0P/AI2O3 samples are presented 
in Fig. 2. Five water peaks evolved through the reduction 
process, corresponding to reduction of different materials 
present in the precursor and to different stages in the re- 
duction process. All of the samples, including the alumina 
support, exhibited a peak at 370 K (with a shoulder at 410- 
450 K) arising from physically adsorbed water. The alumina 
support alone exhibited no significant features above 500 K, 
while the molybdenum phosphate-loaded samples had fea- 
tures corresponding to reduction of the deposited species. 
These features were labeled a, 0, y, and 8 in order of their 
appearance with temperature. The a peak occurred at the 
lowest temperatures, was relatively small, and was present 
in all of the samples. The 8 peak was also present in all of 
the samples, occurred at the highest temperatures, and was 
generally large and diffuse. In the lowest loading materials, 
3.5 and 6.8 wt%, the a and 8 peaks constituted the majority 
of the reduction profile, and are attributed to well-dispersed 
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Fig. 2. Temperature programmed reduction profiles for various. MoPO*/ 
A1 2 0 3 samples reduced to 1223 K at 0.083 Ks" 1 (5 Kmin - " 1 ). 
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Fig. 3. Evolution of MoP during the reduction of MoPOx/A^Os (13 wt% 
MoP). Samples were reduced to intermediate temperatures, quenched in he- 
lium, and characterized by in situ CO uptake before passivation, (o) M0O3; 
(y) M0O2; (+) Mo°; (x) MoP, and (*) A1 2 0 3 . 

species on the alumina surface. The y peak first appeared 
as a small feature in the 6.8 wt% material, and was preva- 
lent at higher loadings along with the 0 peak. The 0 and 
y features slightly overlapped the a and 8 features. But the 
P and y peak temperatures were in the range typically ob- 
served for the reduction of bulk molybdenum compounds, 
and will presently be shown to correspond to X-ray visible 
phase changes. 

The reduction of the 13 wt% M0P/AI2O3 material was 
followed by XRD of quenched and passivated intermedi- 
ates. These spectra are shown in Fig. 3. Comparison of Fig. 2 
with Fig. 3 allows association of the reduction peaks a, 0, 
y, and 8 with phase changes observed by X-ray diffrac- 
tion. The alumina-supported MoP precursor differs from the 
X-ray amorphous precursors for unsupported MoP [1] or 
silica-supported MoP [3] because of the visibility of M0O3. 
The formation of M0O3 is attributed to abstraction of the 
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phosphate component by the alumina surface. The take up 
of phosphate by alumina in standard phosphorus-promoted 
catalysts is well established [24,25]. 

The a peak near 675 K resulted in no new features in 
the XRD spectrum, and was therefore attributed to dehy- 
droxylation of well-dispersed molybdate species. The low 
temperature of this feature is thought to arise because the 
strength of the alumina-molybdate bonds weakens the re- 
maining hydroxyl bonds, which are subsequently easier to 
remove. 

The peak and y peaks correspond to the transforma- 
tion of M0O3 to Mo metal through the intermediate M0O2. 
This behavior is typical of the reduction of bulk M0O3 in hy- 
drogen. This X-ray observable process must occur because 
the phosphorus species are unavailable to the free molybde- 
num species. This demonstrates that the alumina fully binds 
the phosphate at this concentration level (13 wt% MoP). 
At higher surface loadings, the phosphate groups saturated 
the alumina surface and become available to associate with 
the free M0O3. Here, the ft and y peaks resulted in XRD 
changes similar to substoichiometric MoPjcO^ {x < 1) ox- 
ides, where the value of x depends on the relative amount of 
P and Mo. Fortunately for the context of this report, catalysts 
prepared at 1 123 K contained only MoP in their X-ray pat- 
terns (regenerated in the 8 peak at higher temperatures), and 
the catalytic activity was attributed to well-dispersed phases 
not observed by XRD. Therefore, detailed characterization 
of the phase behavior of the samples with loading 23 and 
38 wt% at intermediate temperatures was unnecessary. 

Fig. 3 shows that the <5 peak is accompanied by the dis- 
appearance of molybdenum metal and the appearance of 
molybdenum phosphide above 1123 K. This demonstrates 
that the 8 peak is associated with the freeing of phospho- 
rus within the sample, and its return to the molybdenum 
species. Indirectly, this shows that the 8 peak is associated 
with reduction of alumina-bound phosphate groups, and in- 
dicates that phosphorus has a relatively fast diffusivity in 
molybdenum — because the freed phosphorus can penetrate 
the molybdenum metal completely to form pure MoP. 

Fig. 4 shows the development of the in situ CO uptake 
in intermediates from the reduction of 13 wt% M0P/AI2O3 
superimposed on its reduction profile. The CO uptake cor- 
relates reasonably well with the degree of completion of the 
reduction process, in which the dip at the end is likely due to 
sintering/grain growth of particles of the active phase due to 
the high temperature at that stage. 

Table 2 summarizes the CO uptakes obtained for the 
13 wt% sample at the different stages of reduction marked by 
peaks a, y t and 8. It can be seen that the majority of the 
uptake is generated following the 8 peak. Since the 8 peak 
is also the major reduction peak in the low loading sample 
which does not show an X-ray visible MoP phase, this result 
suggests that the measured CO uptake in all the samples is 
associated with well-dispersed molybdenum phosphide held 
tightly on the alumina surface, rather than the X-ray visible 
phases. 
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Fig. 4. Evolution of in situ CO uptake sites of quenched intermediates dur- 
ing the preparation of 13 wt% M0P/AI2O3. 

Table 2 



Evolution of CO uptake with the reduction peaks a, 0, y, and 8 



Peak 


Temperature 


CO uptake, 


X-ray visible 


event 


interval, K 


umol/g, % 


phase change 


Of 


298-673 


4 


None 




673-763 


15 


M0O3 Mo0 2 


y 


763-863 


10 


M0O2 -*• Mo metal 


s 


863-1173 


71 


Mo metal MoP 
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Fig. 5. Reduction profiles for "blank" alumina samples reduced to 1223 K 
at 0.083 Ks" 1 (5 Kmin- 1 ). (a) 1.16 mmolg- 1 P0 X /A1 2 0 3 ; (b) A1P0 4 ; 
(c) l.iemmolg- 1 MoPO. f /Al 2 0 3 ;and(d) 1.16 mmolg" 1 MoO^/Al 2 0 3 . 

Fig. 5 shows reduction profiles of reference materials 
including commercial AIPO4, the A1 2 0 3 support contain- 
ing only ammonium phosphate, and the Al 2 03 support 
containing only ammonium molybdate. The supported ref- 
erence materials were prepared with procedures (impreg- 
nation, calcination, and peptization) and molar loading 
(1.16 mmolg" 1 A1 2 0 3 ) identical to those used for gen- 
erating the 13 wt% M0P/AI2O3 catalyst. Samples con- 
taining the alumina support exhibited a large dehydration 
peak near 373 K. This low temperature behavior is sim- 
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ilar to that of the alumina blank and the supported cata- 
lysts shown in Fig. 2, and is attributed to the release of 
physically adsorbed water condensed in fine pores. Refer- 
ence samples consisting of only P and AI2O3 had smooth 
baselines between 600 and 900 K, and developed peaks 
above 1 100 K. A yellow-brown deposit formed downstream 
on the interior walls of the quartz reactor tube during the 
high-temperature peaks of the P and AI2O3 samples. The 
deposit is probably elemental phosphorus transported as 
volatile reduced species such as P2, P4, or PH3 [26,27], but 
could be oxidized phosphorus transported as volatile P2O5. 
Vapor-phase transport of phosphorus has been reported ear- 
lier for phosphorus-containing materials undergoing high- 
temperature TPR [19]. The molybdenum-containing mate- 
rials, M0/AI2O3 and M0P/AI2O3, both had identifiable a, 
P, and y peaks. These peaks are respectively attributed 
to recombination of hydroxyl ligands from well-dispersed 
molybdates, M0O3 to M0O2 reduction, and M0O2 to Mo 
metal reduction. The M0/AI2O3 had a tail above 900 K, 
attributed to the reduction of strongly held molybdates. In 
contrast to the reference samples, the M0P/AI2O3 sample 
had a large feature, the 8 peak, between 900 and 1200 K. 
This feature in this temperature range only occurred in the 
presence of both phosphorus and molybdenum. X-ray dif- 
fraction indicated a phase change from Mo to MoP accompa- 
nying the 8 peak, demonstrating influx of phosphorus from a 
source outside the molybdenum grains. The 8 peak in the 
M0P/AI2O3 samples included the reduction of aluminum 
phosphate above 1 100 K, but with the reduction tempera- 
ture of phosphorus decreased because of facilitation of the 
reduction process by molybdenum. 

X-ray diffraction spectra showing the presence of MoP 
in samples of different loading reduced to 1223 K at 
0.0833 Ks^ (5 Kmin"" 1 ) are presented in Fig. 6. MoP was 
the only crystalline product observed, becoming visible at 
6.8 wt%, and increasing in intensity with loading. The XRD 
spectra of the catalysts reduced to 1 123 K at 0.0833 Ks" 1 
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Fig. 6. X-ray diffraction results showing increased amounts of MoP with 
loading for samples reduced to 1223 K at 0.083 Ks -1 (5 Kmin" 1 ). 
(x)MoP,and(*) Al 2 0 3 . 



(5 Kmin -1 ) and held for 2 h were similar to those shown in 
Fig. 6. The lattice constants of the supported MoP catalysts 
are reported in Table 3. The lattice parameters of the MoP in 
the catalysts agree well with the literature values, indicating 
that supported MoP of good purity is formed. 

The effect of final reduction temperature (heating rate 
P = 5 Kmin -1 ) on catalytic performance was evaluated for 
the 13 wt% M0P/AI2O3 material. Of four temperatures in- 
vestigated, successful samples were produced only at 1023 
and 1 123 K. Here, the primary criterion for success was the 
presence of sufficient and stable CO uptake for catalytic test- 
ing. A sample reduced at 873 K for 2 h had insufficient ex 
situ CO uptake (=12 umol g~ l ) to load 70 umol of sites, and 
a sample reduced at 1223 K for 0.5 h exhibited significant 
daily decrease in its CO uptake (noticed because of the aver- 
aging procedure used for catalyst samples) and was therefore 
also not used. Although the sample prepared at 1023 K had 
adequate CO uptake, its XRD spectrum revealed the pres- 
ence of both Mo metal and MoP, while the sample reduced 
at 1 123 K contained only MoP. Because of this, most of the 
supported samples used for the characterization and reactiv- 
ity studies were reduced at 1 123 K. 

Table 3 summarizes crystallite sizes of the molybde- 
num phosphide calculated using the Scherrer equation. 
These values correspond to an average over the low index 
(1000), (1010), and (0001) planes. The sizes of the alumina- 
supported MoP crystallites increased slightly with the load- 
ing amount. More dramatic, however, was their large size 
(about double) compared to the bulk material. This is at- 
tributed to grain growth (Ostwald ripening or sintering) of 
the crystallites at the 200 K higher temperature of 1123 K 
used in the preparation of the supported catalysts compared 
to the bulk material, which was prepared with a maximum 
temperature of 923 K. 

Table 3 also reports site densities of surface metal atoms 
calculated from the crystallite sizes of the X-ray visible MoP 
and the wt% loading of the MoP in the sample. The calcula- 
tions were performed assuming spherical or cubic particles 
having surface site densities equal to 9.71 atoms nm -2 . It is 
understood that the crystallite size as given by the Scherrer 

Table 3 

X-ray diffraction results for M0P/AI2O3 catalysts 



Sample 


Lattice 


Lattice 


Crystallite 


Total metal 




parameter, 


parameter, 


size, 2 


site density, b 




a Q (pm) 


c 0 (pm) 


nm 


pmolg -1 


13%MoP/Al 2 0 3 


321 


318 


39 


43 


23% MoP/Al 2 0 3 


322 


319 


40 


75 


39% M0P/AI2O3 


322 


319 


51 


98 


MoP 


322 


318 


19 


680 


MoP reference 0 


322 


319 







a Average size calculated from line broadening in the (1000), (1010), 
and (0001) directions. 

b Site density = surface metal atoms per gram of sample. Calculated us- 
ing the surface metal areal density of MoP (9.7 x 10 18 surface Mo m~ 2 ), 
the MoP fractional weight loading, and the average MoP crystallite size. 

c JCPDS-ICDD Powder Diffraction File 24-771. 
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Fig. 7. Development of in situ CO uptake as a function of loading. The 

samples were reduced to 1223 K at 5 Kmin" 1 in hydrogen flow at pjg g Effect Qf heating fate Qn the temperature-programmed reduction 
670 M mols »g '(HXXWmin g »). forming 23 wt% MoP/ A1 2 0 3 . 



Table 4 

Results of BET surface area and CO chemisorption measurements on 
M0P/AI2O3 samples 



Sample 


Surface area 


In situ 3 


Ex situ b 


Ex situ b 




fresh, c 


CO uptake, 


CO uptake 


CO uptake 




mV 


umolg^ 1 


fresh, c 


spent, d 








umolg" 1 


umolg" 1 


y-Al 2 0 3 


91 a 


1 


0 




6.8 wt% MoP 


75 


70 


26 


20 


13 wt% MoP 


71 


81 


36 


31 


23 wt% MoP 


60 


76 


28 


29 


39 wt% MoP 


49 


60 


36 


26 



a Samples reduced to 1223 K (0 = 5 Kmin -1 , H 2 flow = 650 umol 
g~ 1 ) and then quenched. 

b Pretreated in H 2 at 723 K for 2 h (/3 = 5 Kmin" 1 , H 2 flow = 650 
umolg" 1 ). 

c Catalysts prepared by reduction to 11 23 K for 2 h (0 = 5 K min~ 1 , H 2 
flow = 650 umol g~ 1 ). 
d Spent catalysts washed in hexane and dried in air. 

equation is only an approximation and ignores the presence 
of a distribution of sizes. The surface site density was cal- 
culated from the crystal structure of MoP [4] assuming ex- 
posure of equal numbers of low index planes. It is assumed 
in these calculations that the surface area to volume ratio of 
the MoP crystallites was proportional to 6/d (spherical or 
cube shaped particles), and that the particles had the den- 
sity of bulk MoP [8] (7.50 gem -3 ). The site density of the 
unsupported MoP was much larger than the experimentally 
determined CO uptake of 18 pmolg" 1 [1]. This reveals that 
the surface sites on bulk MoP are blocked in a ratio near 
38:1, similar to results found for WP [28]. The calculated 
site densities were lower for the supported materials than the 
bulk due to a combination of larger MoP particles sizes and 
lower MoP weight fraction in the mixture. 

Table 4 summarizes measurements of surface area, in situ 
CO uptake, and ex situ CO uptake for MoP/ AI2O3 samples. 
The in situ CO uptakes were collected on the series of sam- 



Table 5 



Apparent first-order activation energies of reduction events in the prepara- 
tion of 23 wt% MoP/Al 2 0 3 



Peak 


Activation energy, 




kJmoP 1 


Of 


152 




157 


y 


130* 


8 


160 



* Nonlinear. Estimated using heating rates of 2 K and 5 K min 1 . See 
text. 



pies whose XRD spectra were presented in Fig. 6 and are 
plotted in Fig. 7. The in situ CO uptakes rose from 1 to 
81 (imol per gram of sample as loading increased from 0 
to 13 wt% MoP and then tapered off to 60 jimol per gram 
as loading was further increased to 39 wt%. The decrease at 
high loadings is likely due to blocking of pores by phosphide 
particles and consequent loss of surface area. Because the ex 
situ CO uptake of the fresh catalyst was used as the basis for 
loading the hydroprocessing reactor, these values were aver- 
aged over three measurements. The ex situ measurements on 
all of the catalyst samples were similar, regardless of load- 
ing, and the results on the fresh and spent catalysts show that 
the catalytic treatment led to only minor, if any, reductions 
in the ex situ CO uptake. Comparison of in situ and ex situ 
CO uptakes reveals that only a portion of the uptake sites are 
recovered after air exposure and subsequent rereduction. 

Heating-rate variation between 0.5 and 5 Kmin" 1 was 
investigated on the sample containing 2.3 mmolg" 1 AI2O3 
(23 wt% M0P/AI2O3). The reduction profiles are displayed 
in Fig. 8, and Table 5 reports the first-order activation en- 
ergies calculated for each type of peak. As heating rate de- 
creased, the relative intensity of the ft peak increased, and 
the y peak resolved into three smaller peaks. Furthermore, 
the y peak maxima did not conform to a linear activation 
energy relationship. This behavior is consistent with satura- 
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tion of phosphorus on the alumina surface between 13 and 
23 wt%, and the overlap of M0O2 -> Mo and MoP^O^ 
MoP* processes in the y peak. The origin of increased inten- 
sity of the ft peak is not clear, but could be related to deeper 
reduction of the crystallites in the slower process. 

5.2. Catalytic activity ofMoP/y-AhOs samples 

The catalysts were evaluated for activity in hydrotreating 
of a simulated distillate containing quinoline (2000 wppm 
N) and dibenzothiophene (3000 wppm S) at 643 K and 
3.1 MPa. Importantly, catalysts were loaded into the reac- 
tor using a constant 70 pmol ex situ CO uptake basis. The 
choice of ex situ CO uptake as the basis for comparison 
assumes catalytic activity to be due to surface sites gener- 
ated through reduction at 723 K, and will be discussed later. 
The use of ex situ uptakes collected of air-exposed samples 
ensures that the sites are regenerated during catalytic pre- 
treatment. Table 4 summarizes CO uptakes obtained after 
pretreatment in hydrogen at 723 K for 2 h. 

Reference materials were also investigated for catalytic 
performance. These standards included the reduced support 
alone, a calcined phosphate blank containing 1.16 mmol 
POjc g" 1 AI2O3, and the phosphate blank exposed to the 
reduction procedure used for the majority of M0P/AI2O3 
samples ( 1 123 K for 2 h). An amount of 3.0 g of each blank 
was loaded, as their CO uptakes were negligible. 

The hydrotreating performances of the samples are re- 
ported in Table 6. The blank reference materials had very 
low HDN and HDS conversions, as expected. For the 
M0P/AI2O3 samples there was only a minor effect of load- 
ing on the HDN and HDS performances, which averaged 56 
and 52%, respectively. The amount of hydrogenated quino- 
line products was stable near 30% (except the phosphorus- 
containing blanks), consistent with equilibration between 
quinoline and 1 ,2,3,4-tetrahydroquinoline. Overall, the re- 
sults indicate that a comparison of activity based on CO 
uptakes is reasonable. The single low HDS conversion of 
44% for the 39 wt% sample is likely due to partial pore 
blockage for this high-loading sample. This would affect 
the larger dibenzothiophene more than the quinoline. The 



Table 6 

Hydroprocessing conversions of M0P/AI2O3 catalysts 3 



Sample 


% HDN 


% HYD 


%HDS 


A1 2 0 3 blank 


2.4 


32 


1.1 


PO x /AI 2 0 3 


3.3 


18 


4.7 


P/AI2O3 


3.8 


19 


1.0 


6.8 wt% MoP 


54 


31 


51 


13 wt% MoP 


52 


33 


57 


13 wt%MoP b 


62 


29 


57 


23 wt% MoP 


58 


28 


53 


39 wt% MoP 


54 


33 


44 


a Reduced to 1 123 K (/J = 5 Kmin" 1 
2h. 

b Reduced to 1023 K (P = 5 Kmin" 1 


, H 2 flow = 
, H 2 flow = 


650 umolg -1 ) for 
650 umolg" 1 ) for 



1.5 h. 



activity of the supported M0P/AI2O3 samples was moder- 
ately high, especially in HDN, and can be compared to that 
of a commercial Ni-Mo-S/AhCh, which at the same hy- 
droprocessing conditions had HDN conversion of 38% and 
HDS conversion of 79% [29]. 

Table 7 reports turnover frequencies (moles converted 
per CO uptake site per second) calculated for quinoline 
HDN and dibenzothiophene HDS for MoP supported on alu- 
mina compared with bulk MoP [1] and MoP on silica [3]. 
The silica-supported material shows somewhat lower activ- 
ity when compared on this basis, and the alumina-supported 
material is similar to the bulk. This suggests that the MoP 
active sites formed on the alumina support are similar to the 
bulk material. Still the differences are small, and again indi- 
cate that CO uptakes are a reasonable method for comparing 
different samples. 

Table 4 also reports the ex situ CO uptake values for the 
spent catalysts. For these measurements the samples were 
removed from the catalytic reactor, washed with hexane, 
and dried in air. Subsequently, they were reduced in hydro- 
gen under the same conditions used for the fresh catalysts. 
The ex situ CO adsorption sites persisted following catalytic 
reaction conditions and exposure to air. This remarkable sta- 
bility in CO uptake values before and after reaction is con- 
sistent with the stability in conversion levels noted for these 
catalysts. 

X-ray diffraction traces for fresh and spent 13 wt% 
M0P/AI2O3 catalyst reduced to 1123 K are presented in 
Fig. 9. MoP and AI2O3 were the only crystalline phases ob- 
served, and the catalyst XRD signature was not significantly 
affected by the catalytic reaction. Similar results were found 
for catalysts of other loadings, whose XRD patterns were 
also stable toward the catalytic test, and which resembled 
the spectra shown in Fig. 6. 

For the 13 wt% M0P/AI2O3 catalyst reduced at 1023 K 
(i.e., 100 K lower than the "standard" temperature of 
1123 K), the X-ray diffraction patterns for the fresh and 
spent catalysts contained features of both Mo and MoP, as 
shown in Fig. 10. Nevertheless, this material was also ac- 
tive and stable toward the catalytic reaction conditions. This 
catalyst was prepared at 1023 K using a 1.5-h soak. Com- 
pared to the intermediate sample quenched at 1023 K, which 
contained only Mo metal in its XRD spectrum (Fig. 3), the 
diffraction trace of the catalyst sample reduced at 1023 K 
contained mostly MoP. Thus, the molybdenum phosphide in 
the catalyst sample appeared during the 1 .5-h soak time en- 



Table 7 

Comparison of hydroprocessing turnover rates of 13% M0P/AI2O3 with 
13% MoP/Si0 2 and bulk MoP 



Sample 


HDN, 


HDS, 




mot quinoline 


mol dibenzothiophene 




mol site" 1 s~ 1 


mol site -1 s -1 


13 wt% MoP/Al 2 0 3 


1.2 x 1CT 3 


8.5 x 10~ 4 


13 wt% MoP/Si0 2 


7.2 x 10^ 4 


3.0 x 10~ 4 


Bulk MoP 


1.5 x 10" 3 


5.2 x 10- 4 
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Degrees / 26 

Fig. 9. X-ray diffraction of fresh and spent 13 wt% MoP/Al 2 0 3 catalyst 
prepared at 1123 K. 




0 1 0 20 30 40 50 60 70 80 90 1 00 
Degrees / 29 

Fig. 10. X-ray diffraction of fresh and spent 13 wt% M0P/AI2O3 catalyst 
prepared at 1023 K. 

tirely by the isothermal progress of the process giving the 8 
peak. 

3.3. Nature of active sites 

The choice of ex situ CO uptake as the basis of compar- 
ison of catalysts deserves commentary. Ordinarily for bulk 
materials, a surface area basis is chosen. However, in sup- 
ported systems the surface area of the active phase and the 
support material combines and the relative surface of the ac- 
tive material is masked. By choosing a surface site titration 
technique (CO uptake) that is known to be roughly propor- 
tional to the total surface area [28], an estimate can be made 
of the amount of active surface available on a supported ma- 
terial. The CO uptake basis then allows comparisons of sup- 
ported materials on the basis of equal active surface. While 
the correlation between CO uptake and catalytic perfor- 
mance is subject to experimental uncertainties, the relative 
correlation between CO uptake and catalytic activity (0. 1 1 



HDN, 0.25 HDS) in this study is substantially better than 
that found between catalytic activity and weight (0.32 HDN, 
0.45 HDS), surface area (0.5 1 HDN, 0.33 HDS), or moles of 
molybdenum (L33 HDN, 1.35 HDS). Here, the relative cor- 
relations represent variabilities, and are defined as the max- 
imum span in values (maximum minus minimum) divided 
by the average value. For a previously undescribed catalyst 
system exhibiting a wide range of loadings and preparation 
conditions, we feel that the consistency of results based on 
CO uptakes is very good. Furthermore, the similarity of the 
turnover frequencies on CO uptake basis with bulk MoP is 
notable. 

The majority of CO uptake in the M0P/AI2O3 system 
was generated by the 8 reduction peak at much higher tem- 
peratures, about 200 K higher, than required for the forma- 
tion of bulk MoP. The 8 peak was observed in all of the 
M0P/AI2O3 samples, even at the lowest loading level of 
3.2 wt%, and is attributed to highly dispersed molybdenum 
and phosphorus on the alumina surface. Our earlier EXAFS 
measurements [2] identified this dispersed phase in reduced 
6.8 wt% M0P/AI2O3 to be clusters of MoP between 0.5 
and 1.2 nm in size. The surface density (Table 1) of MoP 
on the 6.8 wt% sample is 3.8 MoP per nm 2 , and appears 
to be slightly above the saturation limit of molybdenum on 
the alumina in this system (note the presence of traces of 
X-ray visible MoP in Fig. 6, and a slight gamma peak in 
Fig. 2). This saturation limit is slightly less than that of 
M0O3, which forms a saturated monolayer on alumina at 
around 4.5 molecules/nm 2 [30]. Higher surface loadings of 
Mo and P up to 25 MoP/nm 2 (39 wt%) extend well beyond 
monolayer formation on the surface, and the excess molyb- 
denum forms X-ray visible M0O3 until phosphorus satu- 
rates as well, and substoichiometric molybdenum phosphate 
(MoPjcOy, x < 1) species form at higher loadings. These 
X-ray visible excess molybdenum species undergo reduc- 
tion in the temperature range of 650-850 K (at 5 Kmin" 1 ), 
as found in ordinary bulk Mo compounds, but the particles 
are generally very large. The relative sizes of the particles 
(Table 3) are much larger than bulk MoP produced at lower 
temperatures. The sizes of the particles and their dilution by 
alumina lead to much lower values of surface metallic den- 
sities than those found in bulk MoP. The relative similarity 
of turnover frequencies (Table 7) on the basis of CO uptake 
sites in alumina and bulk MoP materials indicates that the 
crystalline phases in the supported catalysts provide only a 
minimal contribution to the active surface in these materials. 

Fig. 1 1 depicts a diagram of the proposed structure of the 
6.8 wt% M0P/AI2O3 and 13 wt% M0P/AI2O3 materials. 
Here, the alumina particles and the associated oxide phases 
are highly simplified for the purpose of clarity. Neverthe- 
less, the salient points of the process are clearly described. 
In the 6.8 wt% loading material, well-dispersed MoPO* re- 
duces directly into MoP particles with diameters on the or- 
der of 1.2 nm. In the 13 wt% MoP material the precursor 
has distributed the molybdenum and phosphorus oxides as 
well dispersed MoPO x , well dispersed aluminum phosphate, 
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Particles 




Fig. 1 1 . Proposed structural model for 7 and 13 wt % M0P/AI2O3 catalysts. 

and crystalline M0O3. During the initial stages of reduction 
(i.e., the £ and y peaks), the crystalline M0O3 reduces di- 
rectly into metal, just as ordinary crystalline M0O3 will do 
in the presence of hydrogen. During the latter stages of re- 
duction (i.e., the 8 peak), the dispersed aluminum phosphate 
is reduced to release phosphorus which is taken up by the 
molybdenum metal, and the well-dispersed MoPO* reduces 
directly into MoP particles with diameters on the order of 
1.2 nm, as was found for the low loading sample. 

3.4. Relationship to previous results in the literature 

Although the nomenclature "M0P/AI2O3" has appeared 
before in the literature [31], the materials studied were in fact 
phosphate-promoted molybdenum sulfides. Such phosphate- 
promoted sulfidic hydrotreating catalysts have been widely 
investigated, as described in an extensive review by Iwamoto 
and Grimblot [13], but are not the same materials investi- 
gated in the present study. In our catalysts, the phosphorus 
is in a reduced state and the active molybdenum phosphide 
is metallic. Aside from the presence of Ni or Co promoters, 
comparison of the preparative conditions reveals that the rel- 
ative amount of P to Mo in the promoted system is small, 
whereas in the present system it has exactly 1 : 1 molar pro- 
portions. For example, optimal activity for thiophene HDS 
in an alumina-supported sulfide catalyst containing 3.0 wt% 
NiO and 15 wt% M0O3 was achieved with a phosphorus 
content of 1 wt% (0.24 mol P: mol Mo) [32]. Furthermore, 
the activation procedure for the promoted system typically 
involved low-temperature (ca. 723 K) reduction in a sulfidic 
atmosphere, whereas in the present study the reduction is 
at high temperature (1 123 K) in pure hydrogen. No defini- 
tive evidence for phosphide formation has been given for the 
promoted system, whereas this study clearly demonstrates 
the presence of hexagonal MoP by XRD. Related work using 
EXAFS shows small crystallites containing Mo-P bonds [2], 

There are also similarities between the phosphide and 
phosphorus-promoted sulfide systems. Both materials ex- 
hibit a strong interaction between phosphate and alumina 
[24,25]. Additionally, the presence of phosphate results in 
moderate reduction in the surface area of the catalysts. It 
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has been found that the interaction of alumina and phosphate 
increases the dispersion of the metals in the Ni-Mo-S sys- 
tem [33]. Molybdena also interacts strongly with alumina, 
suggesting that the phosphate group alters the interaction of 
molybdenum with the surface [33]. 

Because the catalytic activity is roughly independent of 
loading from 6.8 to 39%, while the XRD of MoP increases 
dramatically in intensity, it is reasonable to conclude that 
the majority of active sites are associated with highly dis- 
persed MoP clusters in close contact with AI2O3, as noted 
earlier [2], The nature of phosphorus in M0P/AI2O3 cat- 
alysts reported here is quite different from its oxidized 
state in commercial phosphorus-promoted sulfides. The 
M0P/AI2O3 catalysts reported in this paper are reduced at 
substantially higher temperatures than those employed in 
standard sulfiding, and result in the reduction and transfer 
of phosphorus from alumina to the metallic phase. That the 
temperature of preparation is 200 K higher than for ordinary 
bulk MoP is also remarkable, because it reveals that the alu- 
mina support stabilizes very small MoP particles, protecting 
them from sintering and maintaining their surface sites. 



4. Conclusions 

A series of MoP catalysts supported on /-AI2O3 was 
prepared and tested for hydroprocessing activity under sim- 
ulated industrial conditions. M0P/AI2O3 catalysts were ef- 
fective at HDN and HDS in our tests, averaging 57 and 
52%, respectively. Interestingly, dispersion of the molybde- 
num phosphate precursor on alumina revealed substantially 
different phase behavior than found in bulk MoP and silica 
supported MoP. At loadings between 0 and ~ 6 wt%, both 
the molybdenum and phosphorus species were distributed 
on the alumina support with a strong interaction. At about 
6 wt%, the molybdenum species saturated on the alumina 
support and increases in surface loading led to free M0O3 in 
the oxidized system. Above 13 wt%, the phosphorus also 
saturates on the alumina surface and further increases in 
surface loading produced combinations of M0O3 and sub- 
stoichiometric MoP^O^ (x < 1) species. Regardless of the 
X-ray visible phases, catalytic activity was independent of 
loading and was attributed to well-dispersed MoP formed 
from the reduction of the first 6 wt% of molybdenum phos- 
phate deposited to the surface. 

Higher loadings also resulted in strong interaction be- 
tween alumina and phosphate. Evidence for this was found 
in the presence of M0O3 in the oxidic precursor form of 
the 13 wt% M0P/AI2O3 catalyst, and by new peaks in 
the TPR profile. As reduction proceeded, M0O2 occurred 
at 763 K, molybdenum metal was seen between 863 and 
1073 K, and MoP was formed above 1 173 K. The catalysts 
were active for hydroprocessing of both dibenzothiophene 
and quinoline, and were stable under hydroprocessing con- 
ditions. Analysis suggests that the active sites are associated 
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with small (y 1.2 nm) phosphide crystallites on the alumina 
carrier surface. 
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